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Abstract 

It is shown why the discriminant of a maximal order within a cyclic division algebra must be minimized in 
order to get the densest possible matrix lattices with a prescribed nonvanishing minimum determinant. Using results 
from class field theory a lower bound to the minimum discriminant of a maximal order with a given center and 
index (= the number of Tx/Rx antennas) is derived. Also numerous examples of division algebras achieving our 
bound are given. E.g. we construct a matrix lattice with QAM coefficients that has 2.5 times as many codewords as 
the celebrated Golden code of the same minimum determinant. We describe a general algorithm due to Ivanyos and 
Ronyai for finding maximal orders within a cyclic division algebra and discuss our enhancements to this algorithm. 
We also consider general methods for finding cyclic division algebras of a prescribed index achieving our lower 
\ bound. 

Index Terms 

Cyclic division algebras, dense lattices, discriminants, Hasse invariants, maximal orders, multiple-input multiple- 
^ ' output (MIMO) channels, multiplexing, space-time block codes (STBCs). 

o ' 

I. Overview 

Multiple- antenna wireless communication promises very high data rates, in particular when we have 
^ perfect channel state information (CSI) available at the receiver. In [1] the design criteria for such systems 
O were developed, and further on the evolution of space-time (ST) codes took two directions: trellis codes 
^ ■ and block codes. Our work concentrates on the latter branch. 

£ — ! We are interested in the coherent multiple input-multiple output (MIMO) case. A lattice is a discrete 
O finitely generated free abelian subgroup L of a real or complex finite dimensional vector space V, called 
^ ■ the ambient space. In the space-time setting a natural ambient space is the space M n (C) of complex nxn 
•• . matrices. We only consider full rank lattices that have a basis xi, x-i, . . . , x 2n 2 consisting of matrices that 
.£h ! are linearly independent over the field of real numbers. We can form a 2n 2 x 2n 2 matrix M having rows 
^ I consisting of the real and imaginary parts of all the basis elements. It is well known that the measure, 
■ or hypervolume, m(L) of the fundamental parallelotope of the lattice then equals the absolute value of 
det(M). Alternatively we may use the Gram matrix 

G(L) = MM T = (mr(x l xf)) 1 ^ 2n2 , 

where H indicates the complex conjugate transpose of a matrix. The Gram matrix then has a positive 
determinant equal to m(L) 2 . 

From the pairwise error probability (PEP) point of view [2], the performance of a space-time code 
is dependent on two parameters: diversity gain and coding gain. Diversity gain is the minimum of the 
rank of the difference matrix X — X' taken over all distinct code matrices X, X' E C, also called the 
rank of the code C. When C is full-rank, the coding gain is proportional to the determinant of the matrix 
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(X — X') (X — X') H . The minimum of this determinant taken over all distinct code matrices is called the 
minimum determinant of the code C. If it is bounded away from zero even in the limit as SNR — > oo, the 
ST code is said to have the nonvanishing determinant (NVD) property [3]. For non-zero square matrices, 
being full-rank coincides with being invertible. 

The data rate R in symbols per channel use is given by 

fl=itog| 5 |(|C|), 

where \S\ and \C\ are the sizes of the symbol set and code respectively. This is not to be confused with 
the rate of a code design defined as the ratio of the number of transmitted information symbols to the 
decoding delay (equivalently, block length) of these symbols at the receiver for any given number of 
transmit antennas using any complex signal constellations. If this ratio is equal to the delay, the code is 
said to have full rate. 

The very first STBC for two transmit antennas was the Alamouti code [4] representing multiplication 
in the ring of quaternions. As the quaternions form a division algebra, such matrices must be invertible, 
i.e. the resulting STBC meets the rank criterion. Matrix representations of other division algebras have 
been proposed as STBCs at least in [5]-[14], and (though without explicitly saying so) [15]. The most 
recent work [7]-[15] has concentrated on adding multiplexing gain, i.e. multiple input-multiple output 
(MIMO) applications, and/or combining it with a good minimum determinant. It has been shown in [14] 
that CDA-based square ST codes with the NVD property achieve the diversity-multiplexing gain (D-MG) 
tradeoff introduced in [16]. The codes proposed in this paper all fall into this category and are in that 
sense optimal. Furthermore, algebras with an imaginary quadratic field as a center yield lattices with a 
good minimum determinant, as the corresponding rings of integers have no short non-zero elements. 

Here, yet another design criterion is brought into the playground, namely an explicit criterion for 
maximizing the density of the code. The field of ST coding seems to be lacking a general, precise notion 
for the density in the case of noncommutative structures. In fact, according to our best knowledge the 
theory of orders required for giving this notion has never been considered before in this area. 

Hence, after a cyclic division algebra has been chosen, the next step is to choose a corresponding lattice, 
or what amounts to the same thing, to choose an order within the algebra. Most authors [15], [14] have 
gone with the so-called natural order (see the next section for a definition). One of the points we want 
to emphasize in this article is to use the maximal orders instead. The idea is that one can sometimes use 
several cosets of the natural order without sacrificing anything in terms of the minimum determinant. So 
the study of maximal orders is clearly motivated by an analogy from the theory of error correcting codes: 
why one would use a particular code of a given minimum distance and length, if a larger code with the 
same parameters is available. The standard matrix representation of the natural order results in codes that 
have a so-called threaded layered structure [17]. When a maximal order is used, the code will then also 
extend 'between layers'. However, our simulations suggest that restoring the layered structure somewhat 
by replacing the maximal order with its smartly chosen ideal yields codes with better performance. For 
more details about this see Section IV 111 below. Earlier we have successfully used maximal orders in a 
construction of some 4Tx antenna MISO lattices [5]. 

In some cases the index of the natural order as a sublattice of a maximal order is quite large. E.g. 
in the cases of a family of cyclic algebras suggested in [11] one can theoretically increase the data rate 
by 1.5, 6.5 and 20.5 bits per channel use for 2, 4 and 8 antenna codes, respectively. We do emphasize 
that such increments of data rates are only theoretical in nature. This is because one is compelled to use 
relatively large subsets of the infinite lattice before the full density advantage of the maximal order is 
attained. Also the lattice of a fully multiplexing 8Tx+8Rx antenna MIMO code has dimension 128. The 
nearest vector problem in such high-dimensional lattices is used in some cryptographic applications, so it 
is safe to say that ML-decoding of such lattices will have prohibitive complexity. These numbers, however, 
motivated us to look for methods of locating maximal orders. A general purpose algorithm for this task 
has been developed by Ivanyos and Ronyai [18]. A commercially available version of their algorithm is 
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implemented by W. van de Graaf as part of the computer algebra system MAGMA [19]. It turned out 
that this general purpose algorithm was not able to handle the algebras of index eight. To deal with these 
special cases we developed some enhancements to their algorithm. 

Given that maximal orders provide the best codes in terms of minimum determinant vs. average power 
we are left with the question: Which division algebra should we use? To continue the analogy from the 
theory of error-correcting codes we want to find the codes with the highest possible density. That is, with 
the smallest fundamental parallelotope. To that end we need a suitable tool for parameterizing the cyclic 
division algebras with a given center and index. Luckily, relatively deep results from class field theory 
provide us with the necessary tool of Hasse invariants. The measure of a fundamental parallelotope of 
a maximal order (that will later on be referred to as the discriminant of the division algebra) can be 
expressed in terms of Hasse invariants [20]. With these results at hand we then derive a lower bound to 
the discriminant. While the proof of the lower bound is not constructive per se, it does show that our 
lower bound is achievable. In the latter parts of this article we describe some techniques for constructing 
division algebras with a minimal discriminant. 

It is worth mentioning that in [21] the authors have made a similar approach in the reduced case of 
commutative number fields. 

While our interest in these problems is mostly theoretical, some of the densest lattices we have found 
also perform well in computer simulations. Our construction of the densest 2x2 matrix lattice improves 
upon the deservedly celebrated Golden code in block error rates by about 0.9 dB at data rates from 5 
to 6 bpcu. The performance of both the rival codes can be further improved by coset optimization and 
this also cuts down the gap to about 0.3 dB. Observe that at the data rate of 4 bpcu we have a tie. This 
is easily explained by the fact that for codes of that size there is a particularly attractive choice for the 
coset of the Golden code. Another point worth keeping in mind is that the somewhat irregular geometry 
of our lattice more or less necessitates the use of a code book as opposed to a simple combination of 
Gray coding and PAM. However, this also holds for the Golden code, when we do any coset optimization. 
Thus we might conclude that our work shows that not using a codebook costs about 1 dB. 

The paper is organized as follows. In Section [III various algebraic notions related to cyclic algebras, 
Brauer groups, orders, discriminants, and localizations are introduced and demonstrated by examples. 
Furthermore, it is shown that maximizing the density of the code, i.e. minimizing the fundamental 
parallelotope is equivalent to minimizing the discriminant. This leads us to Section [Till where we derive an 
achievable lower bound for the discriminant. In Section [TV] we propose a general algorithm due to Ivanyos 
and Ronyai [18] for finding maximal orders. Unfortunately, when we were trying to use the MAGMA 
implementation of this algorithm for finding maximal orders of certain cyclic division algebra of index 
no more than 8, the memory of a typical modern PC turned out to be insufficient. Hence, also some 
enhancements to their algorithm in this special case are discussed here. The Perfect codes are analyzed 
in Section [V] in terms of Hasse invariants and discriminants. We show that the natural orders (i.e. the 
orders the authors have used in [10]) of the related algebras are maximal in the cases of #Tx = 2 and 
#Tx = 3, but can be enlarged in the cases of #Tx = 4 and #Tx = 6. In Section [VI] we construct 
division algebras with a minimal discriminant. The case of a unit non-norm element is separated from 
the general construction. Finally in Section IVIIl the theory is brought into practice by giving an explicit 
code construction that outperforms or ties with the Golden code. Simulation results are provided to back 
up this claim. 

II. Cyclic algebras, Brauer groups, orders, and discriminants 

We refer the interested reader to [22] and [7] for a detailed exposition of the theory of simple algebras, 
cyclic algebras, their matrix representations and their use in ST-coding. We only recall the basic definitions 
and notations here. In the following, we consider number field extensions E/F, where F denotes the 
base field and F* (resp. E*) denotes the set of the non-zero elements of F (resp. E). In the interesting 
cases F is an imaginary quadratic field, either Q(i) or Q(\/— 3). We assume that E/F is a cyclic field 
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extension of degree n with Galois group Gal(F/F) = (a). Let A = (F/F, (7,7) be the corresponding 
cyclic algebra of degree n (n is also called the index of A), that is 

A = E © uE © u 2 E © • • • © u n - l E, 

as a (right) vector space over E. Here u G .4 is an auxiliary generating element subject to the relations 

xu = ua(x) for all x G E and u n = 7 G F*. An element a = x + -ua;i + • ■ • + ■u n_1 a; n _ 1 G .4 has the 
following representation as a matrix A = 

/ x 7cr(x„_i) 7(T 2 (a; n _2) ••• 7^ 1 ( a; i)\ 
Xi cr(x ) 7cr 2 (a; n _i) 'ya n ~ 1 (x 2 
x 2 o-(xi) o- 2 (x ) 7<T n_1 (a;3 

o-(x n _ 2 ) (T 2 (a;n-3) ••• o-™" 1 ^) / 

We refer to this as the standard matrix representation of A. Observe that some variations are possible 
here. E.g. one may move the coefficients 7 from the upper triangle to the lower triangle by conjugating 
this matrix with a suitable diagonal matrix. Similarly one may arrange to have the first row to contain the 
"pure" coefficients x , . . . , x n ^\. Such changes do not affect the minimum determinant nor the density of 
the resulting lattices. 

If we denote the basis of E over Of by {1, ei, e„_i}, then the elements x i: % — 0, n — 1 in the 
above matrix take the form xi = Ylk=o fk&k, where fk G Of for all k = 0, n — 1. Hence n complex 
symbols are transmitted per channel use, i.e. the design has rate n. In literature this is often referred to 
as having a. full rate. 

Definition 2.1: The determinant (resp. trace) of the matrix A above is called the reduced norm (resp. 
reduced trace) of the element a G A and is denoted by nr(a) (resp. tr(a)). 

Remark 2.1: The connection with the usual norm map Na/f(o>) (resp. trace map Ta/f(o>)) and the 
reduced norm nr(a) (resp. reduced trace tr(a)) of an element a G A is N^/f(o) = [nr{a)) n (resp. 
T A / F (a) = ntr(a)), where n is the degree of E/F. 

Definition 2.2: An algebra A is called simple if it has no nontrivial ideals. An F-algebra A is central 
if its center Z(A) = {a G A \ aa' = a'a Va' G A} = F. 

Definition 2.3: Let S denote an arbitrary ring with identity. The Jacobson radical of the ring S is the 
set Rad(S) = 

{x G S I xM = for all simple left ^-modules M}. 

Rad(S') is a two-sided ideal in S containing every nilpotent (i.e. for which X k = for some k G Z+) 
one-sided ideal X of S*. Also, Rad(S') can be characterized as the intersection of the maximal left ideals 
in S. If S is a finite dimensional algebra over a field or, more generally, left or right Artinian then Rad(S') 
is the maximal nilpotent ideal in S. 

A division algebra may be represented as a cyclic algebra in many ways as demonstrated by the 
following example. 

Example 2.1: The division algebra QA used in [3] to construct the Golden code is gotten as a cyclic 
algebra with F = Q(i), E = Q(i, y/5), 7 = 2, when the F-automorphism a is determined by a(y/E) = 
— \fb. We also note that in addition to this representation QA can be given another construction as a 
cyclic algebra. As now u 2 = i we immediately see that F(u) is a subfield of QA that is isomorphic to the 
eighth cyclotomic field E' = Q((), where £ = (1 + i)/y/2. The relation u\^E = —\^Eu read differently 
means that we can view u as the complex number ( and \/5 as the auxiliary generator, call it u' = V5. 
We thus see that the cyclic algebra 

E' © u'E' — (E'/F, a', 7') 

is isomorphic to the Golden algebra. Here a' is the F-automorphism of E' determined by ^ 1 — > — and 

i = u' 2 = 5. 
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Any cyclic algebra is a central simple F-algebra (cf. Definition I2.2I ). Two central simple F-algebras 
A and B are said to be similar, if there exist integers m an n such that the matrix algebras Ai n (A) and 
A4 m (B) are isomorphic F-algebras. Wedderburn's structure theorem [22, Theorem, p. 171] tells us that 
any central simple algebra is a matrix algebra over a central simple division algebra, and it easily follows 
that within any similarity class there is a unique division algebra. Similarity classes of central simple 
algebras form a group (under tensor product over F), called the Brauer group Br(F) of the field F. If 
F' is an extension field of F, and A is a central simple F-algebra, then the tensor product A' = A ®p F' 
is a central simple F'-algebra. We refer to this algebra as the algebra gotten from A by extending the 
scalar s to F' . 

The next proposition due to A. A. Albert [23, Theorem 11.12, p. 184] tells us when a cyclic algebra 
is a division algebra. 

Proposition 2.1 (Norm condition): The cyclic algebra A = (E/F,a, 7) of degree n is a division algebra 
if and only if the smallest factor t E Z + of n such that 7* is the norm of some element of E* is n. 

Due to the above proposition, the element 7 is often referred to as the non-norm element. 

Let F be an algebraic number field that is finite dimensional over Q. Denote its ring of integers by Op. 
If P is a prime ideal of Op, we denote the P-adic completion of F by Fp. The division algebras over F P 
are easy to describe. They are all gotten as cyclic algebras of the form A(n, r) = (E / F P , a,n r ), where 
E is the unique unramified extension of Fp of degree n, a is the Frobenius automorphism, and n is the 
prime element of Fp. The quantity r/n is called the Hasse invariant of this algebra and n is referred to 
as the local index. It immediately follows from Proposition 12.11 that A(n, r) is a division algebra, if and 
only if (r, n) = 1. For a description of the theory of Hasse invariants we refer the reader to [20, p. 266] 
or [24]. 

We are now ready to present some of the basic definitions and results from the theory of maximal 
orders. The general theory of maximal orders can be found in [20]. 

Let R denote a Noetherian integral domain with a quotient field F, and let A be a finite dimensional 
F-algebra. 

Definition 2.4: An R-order in the F-algebra A is a subring A of A, having the same identity element 
as A, and such that A is a finitely generated module over R and generates A as a linear space over F. 
An order A is called maximal, if it is not properly contained in any other _R-order. 

Let us illustrate the above definition by concrete examples. 

Example 2.2: (a) Orders always exist: If M is a full i?-lattice in A, i.e. FM = A, then the left order 
of M defined as 0\{M) = {x E A \ xM C M} is an i?-order in A. The right order is defined in an 
analogous way. 

(b) If R is the ring of integers Op of the number field F, then the ring of integers Op of the extension 
field E is the unique maximal order in E. For example, in the case of the cyclotomic field E = Q(C), 
where ( = exp(2iri/k) is a primitive root of order k the maximal order is Op = Z[£]. 

The next proposition (see [26, proof of Theorem 3.2]) is useful when computing left orders in Section 
H3 

Proposition 2.2: Let A be a simple algebra over F and M a finitely generated Op -module such that 
FM = A. Then there exists an element s E O f \ {0} such that s ■ 1 E M. Moreover, O t (M) = {b E 
s~ l M \bM <M} < s^M. 

For the purposes of constructing MIMO lattices the reason for concentrating on orders is summarized 
in the following proposition (e.g. [20, Theorem 10.1, p. 125]). We simply rephrase it here in the language 
of MIMO-lattices. We often (admittedly somewhat inaccurately) identify an order (or its subsets) with its 
standard matrix representation. 

Proposition 2.3: Let A be an order in a cyclic division algebra (E/F, a, 7). Then for any non-zero 
element a E A its reduced norm nr(a) is a non-zero element of the ring of integers Op of the center F. 
In particular, if F is an imaginary quadratic number field, then the minimum determinant of the lattice A 
is equal to one. 
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Example 2.3: In any cyclic algebra we can always choose the element 7 e F* to be an algebraic 
integer. We immediately see that the C^-module 

A = e> £ © m0 e © • • • © n^ 1 ^, 

where O e is the ring of integers, is an O^-order in the cyclic algebra (E/F, <r, 7). We refer to this 
Oi? -order as the natural order. It will also serve as a starting point when searching for maximal orders. 

We want the reader to note that in any central simple algebra a maximal Z-order is a maximal C^-order 
as well. Note also that if 7 is not an algebraic integer, then A fails to be closed under multiplication. This 
may adversely affect the minimum determinant of the resulting matrix lattice, as elements not belonging 
to an order may have non-integral (and hence small) norms. 

We remark that the term 'natural order' is somewhat misleading. While it is the first order that comes 
to mind, there is nothing canonical about it. Indeed, distinct realizations of a given division algebra as a 
cyclic algebra often lead to different natural orders. E.g. constructing the algebra of rational Hamiltonian 
quaternions from the cyclic extension Q(\/— 3)/Q as opposed to the more common Q(i)/Q leads to a 
different natural order. The interested reader may verify this as an exercise by starting with the observation 
that the Hamiltonian quaternion i + j + k may be used as a square root of —3. 

Definition 2.5: Let m = dim F A. The discriminant of the i?-order A is the ideal d(A/R) in R generated 
by the set 

{det(tr(xiXj))™ j=1 | (x u ...,x m ) G A m }. 
In the interesting cases of F = Q(i) (resp. F = Q(\/— 3)) the ring R = Z[i] (resp. R = Z[u], 
to = (—1 + y/—3)/2) is a Euclidean domain, so in these cases (as well as in the case R = Z) it makes 
sense to speak of the discriminant as an element of R rather than as an ideal. We simply pick a generator 
of the discriminant ideal, and call it the discriminant. Equivalently we can compute the discriminant as 

d(A/R) = &et(tr(xiXj))™ j=1 , 

where {x 1 , . . . ,x m } is any i?-basis of A. It is readily seen that whenever A C T are two i?-orders, then 
d(T) is a factor of d(A). The index [r : A] is related to discriminants by the following lemma. 
Lemma 2.4: 

[R : d(A)R] = [r : A] 2 [R : d(T)R] 
Proof: [20, p.66] ■ 
It turns out (cf. [20, Theorem 25.3]) that all the maximal orders of a division algebra share the same 
discriminant that we will refer to as the discriminant of the division algebra. In this sense a maximal 
order has the smallest possible discriminant among all orders within a given division algebra, as all the 
orders are contained in the maximal one. 

For an easy reference we also note the following basic formula for the discriminant of certain cyclotomic 
fields. 

Proposition 2.5: Let Q = exp(27ri/2 £ ) be a complex primitive root of unity of order 2 £ , where £ > 2 
is an integer. Then n = [Q((e) ■ Q(«)] = 2^ 2 and 

d(z[0]/z[*]) = (i + *) 2ri( '- 2) - 

The definition of the discriminant closely resembles that of the Gram matrix of a lattice, so the following 
results are unsurprising and probably well known. We include them for lack of a suitable reference. 

Lemma 2.6: Assume that F is an imaginary quadratic number field and that 1 and 9 form a Z-basis of 
its ring of integers R. Assume further that the order A is a free i?-module (an assumption automatically 
satisfied, when R is a principal ideal domain). Then the measure of the fundamental parallelotope equals 

m(A) = \%ef\d{A/R)\. 

Proof: Let A = (a^) be an n x n complex matrix. We flatten it out into a 2 x 2n 2 matrix L(A) 
by first forming a vector of length n 2 out of the entries (e.g. row by row) and then replacing a complex 
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number z by a diagonal two by two matrix with entries z and z* (= the usual complex conjugate of z). 
If A and B are two square matrices with n rows we can easily verify the identities 

and 

jT\T _ (tr(AB) 



W(B J ) J = ^ ^ ' (r(j4B) ,j. (2) 

Next let £> = {xi, X2, . . . , x n 2} be an i?-basis for A. We form the 2n 2 x 2n 2 matrix L(B) by stacking 
the matrices L(xi) on top of each other. Similarly we get R(B) by using the matrices L(xf) T as 'column 
blocks'. Then by © the matrix M = L(B)R(B) consists of two by two blocks of the form 



L(xi)L(xT) T 



tr(xiXj) 
tr(xiXjY 

Clearly det R(B) = ±det L(B), and det M = \d(A/R)\ 2 , so we get 

|d(A/i2)| = |detL(£)|. 

Next we turn our attention to the Gram matrix. By our assumptions the set B U 9B is a Z-basis for A. 
Let us denote 

D= 1 1 
\9 9* 

From the identities $t(xy*) = (xy* + x*y)/2 and 

fx 0\ _ ( x x* 
\0 x*J ~ \9x 9*x* 

together with ([T]) it follows that for any two n x n matrices A and B we have 

1 (r>T(A\\(r>T(-R\\ H - f^tr(AB H )) §l(tr(A(9B) H ) 

2 {UL{A)) t VL ^>) ~ {jft(tr(9AB H )) $t(tr(9A(9B) H )) 

Therefore, if we denote by the 2n 2 x 2n 2 matrix having n 2 copies of D along the diagonal and 
zeros elsewhere, we get the following formula for the Gram matrix 

G(A) = i (DWL(iB)) (D^L(B)) H . 



Thus, 

m(A) =detG(A) 1/2 = |detL(fi)| 



-detD 
2 



Our claim now follows from all these computations and the fact that (det D)/2 = [9* — 9)/2 = —^s9. ■ 
In the respective cases F = Q(z) and F = Q(-\/— 3) we have 9 = i and 9 = (— 1 + 3)/2 respectively, 

so we immediately get the following two corollaries. 

Corollary 2.7: Let F = Q(i),i? = Z[i], and assume that A C (E/F,a,j) is an i?-order. Then the 

measure of the fundamental parallelotope equals 

m(A) = \d(A/Z\i])\- 

Example 2.4: When we scale the Golden code [3] to have a unit minimum determinant, all the 8 
elements of its Z-basis will have length 5 1 / 4 and the measure of the fundamental parallelotope is thus 
25. In view of all of the above this is also a consequence of the fact that the Z[z] -discriminant of the 
natural order of the Golden algebra is equal to 25. As was observed in [25] the natural order happens to 
be maximal in this case, so the Golden code cannot be improved upon by enlarging the order within QA. 
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Corollary 2.8: Let u = (-l + v /3 3)/2, F = Q(u), R = Z[uj], and assume that A C (E/F, a, 7) is an 
i?-order. Then the measure of the fundamental parallelotope equals 

m (A) = (v / 3/2)"V(A/ZM)|. 
The upshot is that in both cases maximizing the density of the code, i.e. minimizing the fundamental 
parallelotope, is equivalent to minimizing the discriminant. Thus, in order to get the densest MIMO- 
codes we need to look for division algebras that have a maximal order with as small a discriminant as 
possible. 

For an easy reference we also include the following result that is a relatively easy consequence of the 
definitions. 

Lemma 2.9: Let E/F be as above, assume that 7 is an algebraic integer of F, and let A be the 
natural order of Example 12.31 If d(E / F) is the (^-discriminant of Oe (often referred to as the relative 
discriminant of the extension E/F), then 

d(A/0 F ) = d(£/F) n 7 n(n-1) - 

Proof: In the expansion 

A = Oe © uO E © • • • © u^Oe 

we see that u % Oe and v?Oe are orthogonal to each other with respect to the bilinear form given by the 
reduced trace except in the cases where i + j = (mod n). Assume that i + j is divisible by n for some 
i,j in the range < i, j < n, and that X\, . . . ,x n are elements of Oe- Then the multiplication rules of 
the cyclic algebra imply that 

det(tr(u % XkU^Xi))^ i=1 = ±det(u l+ Hr(xkXf))1 l=1 = ±7 e det(tr(xkXt))^ e=1 , 

where the exponent e is equal to zero or n according to whether i + j equals zero or n. The former case 
occurs only once and the latter case occurs exactly n — 1 times. The claimed formula then follows. ■ 
Example 2.5: We use the notation from Proposition 12.51 In [11] Kiran and Rajan have shown that the 
family of cyclic algebras At = (Q(Q)/Q(i),cr(Q) = 0,2 + with £ > 3, consists entirely of division 
algebras. Let A nat; i be the natural order of the algebra At. We may now conclude from Lemma 12.91 
Proposition 12.5 L and Corollary 12.71 that 

d{Kt,nat/m) = (2 + i) n{n - l \l + lf n ^~ 2 \ 

and that 

m(A nat/ ) 2 = 2 2 ^- 2 ) 5 ^- 1 >. 

For instance, in the 2 antenna case £ = 3,n = 2, we have m(A natt t) = 80, and thus the Golden code is 
denser than the corresponding lattice A3 of the same minimum determinant. However, the natural order 
of A3 is not maximal and we will return to this example later on. 

In Section [IV] some facts from the local theory of orders are required. For the basic properties of 
localization the reader can turn to [22, Chapter 7] or [20, Chapters 1, 2]. For the proofs for the rest of 
this section, see [18] and [26]. 

If R is a Dedekind domain with a quotient field F, and P is a prime ideal in R, then the ring of 
quotients Rp = (R/P)~ 1 R C F is a discrete valuation ring. For the i?-lattices M in A the localization 
at P is defined as M P = R P M C A M P is an i? P -lattice. Moreover, if M is a full (cf. Example |X2l 
i?-lattice in A, then M P is a full -Rp-lattice in A. To be more specific, let us define the ring Z p . 

Definition 2.6: For a rational prime p let Z p denote the ring 

Z P = G Q I r, s G Z, gcd(p, s) = 1}. 

Z p is a discrete valuation ring with the unique maximal ideal pZ p . If A is a Z-order we use the notation 
A p = Z P A. 
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We remark that one should not confuse the localization R P with the ring of integers Pp of the P-adic 
completion. We use the caret to indicate a complete structure. This is somewhat non-standard in the case 
of Z p that is nearly universally used to denote the complete ring of p-adic integers. We use Z p for the 
complete ring. 

The next statement illustrates a simple but useful connection between the orders A and A p . 

Proposition 2.10: Let A be a Z-order in A. The map $ : x i— ► x+pA p , x £ A induces an isomorphism 
of the rings A/pA = A p /pA p . 

Proposition 2.11: Let P be a prime ideal of the ring R. The residue class ring A = A/PA is an 
algebra with identity element over the residue class field R = R/P and dimpA = dim-^A. If : A — > A 
is the canonical epimorphism, then PA C Rad(A) = _1 Rad(A) and <p induces a ring isomorphism 
A/ Rad(A) = A/ Rad(A). As a consequence, a left (or right) ideal X of A is contained in Rad(A) if and 
only if there exists a positive integer t such that X* C PA. 

The following facts establish some practical connections between the local and global properties of 
orders. 

Proposition 2.12: Let A be a simple algebra over F. Let P be a prime ideal of R, and T be an P-order 
in A. Then 

(i) Tp is an Pp-order in A. 

(ii) T is a maximal P-order in A if and only if Tp is a maximal Pp-order in A for every prime ideal 
P of P. 

(iii) d(T/R) P = d{T P /R P ). 

Proposition 2.13: Let P be a prime ideal of P and T be an P-order such that Tp is not a maximal 
Pp-order. Then there exists an ideal X > Pr of T for which 0;(X) > T. 

Extremal orders and especially Proposition 12.151 below play a key role in the method for constructing 
maximal orders. 

Definition 2.7: We say that T P radically contains Ap if and only if Ap C T P and Rad(Ap) C Rad(Tp). 
The orders maximal with respect to this partial ordering are called extremal. Maximal orders are obviously 
extremal. 

Proposition 2.14: An Pp-order Ap is extremal if and only if Ap = C^(Rad(Ap)). 

Proposition 2.15: Let Ap C Tp be Pp-orders in A. Suppose that Ap is extremal and that Tp is minimal 
among the Pp-orders properly containing Ap. Then there exists an ideal J of Ap minimal among those 
containing Rad(Ap) such that Oi(J) D T P . 

III. Discriminant bound 

Again let F be an algebraic number field that is finite dimensional over Q, Of its ring of integers, P 
a prime ideal of Op and F P the completion. In what follows we discuss the size of ideals of Of- By this 
we mean that ideals are ordered by the absolute values of their norms to Q, so e.g. in the case Of = Z[z] 
we say that the prime ideal generated by 2 + i is smaller than the prime ideal generated by 3 as they have 
norms 5 and 9, respectively. 

The following relatively deep result from class field theory is the key for deriving the discriminant 
bound. Assume that the field F is totally complex. Then we have the fundamental exact sequence of 
Brauer groups (see e.g. [20] or [24]) 

— ► Br(P) — ► © Br(Pp) — > Q/Z — ► 0. (3) 

Here the first nontrivial map is gotten by mapping the similarity class of a central division P-algebra V 
to a vector consisting of the similarity classes of all the simple algebras V P gotten from V by extending 
the scalars from F to F P , where P ranges over all the prime ideals of Op. Observe that V P is not 
necessarily a division algebra, but by Wedderburn's theorem [22, p. 203] it can be written in the form 

Vp = M Kp (A P ), 
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where Ap is a division algebra with a center P P , and Kp is a natural number called the /oca/ capacity. The 
second nontrivial map of the fundamental exact sequence is then simply the sum of the Hasse invariants 
of the division algebras Ap representing elements of the Brauer groups Bt(P p ). 

This exact sequence tacitly contains the piece of information that for all but finitely many primes P 
the resulting algebra Vp is actually in the trivial similarity class of Fp-algebras. In other words Dp is 
isomorphic to a matrix algebra over Pp. More importantly, the sequence tells us that the sum of the 
nontrivial Hasse invariants of any central division algebras must be an integer. Furthermore, this is the 
only constraint for the Hasse invariants, i.e. any combination of Hasse invariants (a /nip) such that only 
finitely many of them are non-zero, and that they sum up to an integer, is realized as a collection of the 
Hasse invariants of some central division algebra V over F. 

Let us now suppose that with a given number field F we would like to produce a division algebra A of 
a given index n, having F as its center and the smallest possible discriminant. We proceed to show that 
while we cannot give an explicit description of the algebra A in all the cases, we can derive an explicit 
formula for its discriminant. 

Theorem 3.1: Assume that the field F is totally complex and that Pi, . . . ,P n are some prime ideals of 
Of- Assume further that a sequence of rational numbers ai/m Pl , . . . , a n /mp n satisfies 



n 



E 



mp. 

t=i 1 



= (mod 1), 



1 < Ct>i < m Pv an d ( a i, m Pi) = 1- 

Then there exists a central division F-algebra A that has local indices m Pi and the least common 
multiple (LCM) of the numbers {nip^ as an index. 

If A is a maximal Op-order in A, then the discriminant of A is 

d(A/o F ) = np; 1 ,m *. 

i=l 

Proof: By exactness of the sequence (|3J we know that there exists a central division algebra A over 
F which has local indices m Pi . From [20, Theorem 32.19] we know that ^[A : F] = LCM{m Pi }. By 
[20, Theorem 32.1] the discriminant then equals 

/ n 

d(A/R)= \Y[P t {mp ^ 1)fip A , (4) 

where Kp t is the local capacity. 

A simple calculation of dimensions shows that 



Kp 



nip 

Substituting this into © we get the claim. ■ 
At this point it is clear that the discriminant d(A) of a division algebra only depends on its local indices 

m Pi . 

Now we have an optimization problem to solve. Given the center F and an integer n we should decide 
how to choose the local indices and the Hasse invariants so that the LCM of the local indices is n, the 
sum of the Hasse invariants is an integer, and that the resulting discriminant is as small as possible. We 
immediately observe that at least two of the Hasse invariants must be non-integral. 

Observe that the exponent d(P) of the prime ideal P in the discriminant formula 

d(P) = (nip - l)tii^l = n 2 (l - — 
nip \ nip 
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As for the nontrivial Hasse invariants n > m P > 2, we see that n 2 /2 < d(P) < n(n — 1). Therefore 
the nontrivial exponents are roughly of the same size. E.g. when n — 6, d(P) will be either 18, 24 or 
30 according to whether m P is 2, 3 or 6. Not surprisingly, it turns out that the optimal choice is to have 
only two non-zero Hasse invariants and to associate these with the two smallest prime ideals of Op- 
Theorem 3.2 (Main Theorem): Assume that F is a totally complex number field, and that P\ and P2 
are the two smallest prime ideals in Of. Then the smallest possible discriminant of all central division 
algebras over F of index n is 

(P 1 P 2 ) n(n " 1) . 

Proof: By Theorem 13. II the division algebra with Hasse invariants 1/n and (n— l)/n at the primes 
Pi and P 2 has the prescribed discriminant, so we only need to show that this is the smallest possible 
value. 

By the above discussion it is clear that in order to minimize the discriminant one cannot have more than 
three nontrivial Hasse invariants. This is because for prime ideals Pi, P 2 , P3, Pa (listed from the smallest 
to the largest) we always have 

pd(Pi) pd(P 2 ) pd(Pa) pd(P 4 ) > (p i p 2 )«(n-l) 

as the exponents d(Pi) > n 2 /2 irrespective of the values of the Hasse invariants. A possibility is that 
some combination of three Hasse invariants might yield a smaller discriminant. Let us study this in detail. 

If one of the local indices, say mp 1 , has only a single prime factor, say p, then we can add this Hasse 
invariant together with one of the other two, as long as we are careful to choose the one, say rrip 2 , whose 
denominator is divisible by a smaller power of p. In this addition process the least common multiple of 
the denominators remains the same, so the new set of only two nontrivial Hasse invariants corresponds 
to a division algebra of the same index. This is because in the sum of the Hasse invariants 

ai/mpj + a 2 /mp 2 = a! /m' P , (mod 1) 

the new local index m' pi is gotten from the old local index rrip 2 by multiplying it with a (possibly the 
zeroth) power of p. Let P' be smaller of the two ideals Pi and P 2 . As d(Pi) + d(P 2 ) > n(n — 1) > d'(P'), 
where d'(P') is the exponent corresponding to the local index m P i, this new division algebra (with 
nontrivial Hasse invariants associated with primes P' and P 3 only) will have a smaller discriminant. 

The remaining case is that all the three local indices have at least two distinct prime factors. In this 
case all the three Hasse invariants have numerators > 6. As then d(Pi) + d(P 2 ) + d(P 5 ) > 2n[n — 1), we 
see that the discriminant of the division algebra with these Hasse invariants also exceeds the stated lower 
bound. ■ 

We remark that in the most interesting (for MIMO) cases n = 2 and n — 3, the proof of Theorem 13 .21 is 
more or less an immediate corollary of Theorem 13.11 We also remark that the division algebra achieving 
our bound is by no means unique. E.g. any pair of Hasse invariants a/n, (n — a)/n, where < a < n, 
and (a,n) = 1, leads to a division algebra with the same discriminant. 

The smallest primes of the ring Z[i] are 1 + i and 2 ± i. They have norms 2 and 5 respectively. The 
smallest primes of the ring Z[cu] are -J— 3 and 2 with respective norms 3 and 4. Together with Corollaries 
12.71 and 12.81 we have arrived at the following bounds. 

Corollary 3.3 (Discriminant bound): Let A be an order of a central division algebra of index n over 
the field Q(z). Then the measure of a fundamental parallelotope of the corresponding lattice 

m(A) > io™^- 1 )/ 2 . 

Corollary 3.4 (Discriminant bound): Let A be an order of a central division algebra of index n over the 
field Q(u>), uj = (—1 + a/— 3)/2. Then the measure of a fundamental parallelotope of the corresponding 
lattice 

m(A) > (V?>/2) W2 l2 n{n - l V 2 . 

The Golden algebra reviewed in Example 12.11 has its nontrivial Hasse invariants corresponding to the 
primes 2 + i and 2 — i and hence cannot be an algebra achieving the bound of Theorem 13.21 A clue for 
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finding the optimal division algebra is hidden in the alternative description of the Golden algebra given 
in Example 12.11 It turns out that in the case F = Q(i), E = Q(C) instead of using 7' = 5 as in the case 
of the Golden algebra we can use its prime factor 7 = 2 + i. 

Proposition 3.5: The maximal orders of the cyclic division algebra A3 = (Q(£)/Q(z), cr, 2 + i) of 
Example 12.51 achieve the bound of Theorer riX2l 

Proof: The algebra A3 is generated as a Q(i)-algebra by the elements £ and u subject to the relations 
( 2 = i, u 2 = 2 + i, and u( = —(u. The natural order Z[£] © uZ[(] is not maximal. Let us use the matrix 
representation of A3 as 2 x 2 matrices with entries in Q(C)> so elements of Q(z) are mapped to scalar 
matrices and ( is mapped to a diagonal matrix with diagonal elements £ and We observe that the 
matrix 

If 2i-(l-i) v / 2 (2 + i)(2i- (l + i)y/2)\ 
w ~ 4 V(l+i)(l + V2 + i) 2i + (l-i)V2 / 
is an element of A3. Straightforward calculations show that w satisfies the equations 

w 2 = — i + iw and w( = — 1 + £ 3 — C w - 

From these relations it is obvious that the free Z[£] -module with basis elements 1 and w is an order A. 
Another straightforward computation shows that d(A/Z[i]) = — 8 + 6i = (1 + i) 2 (2 + i) 2 . As this is the 
bound of Theorem 13.21 we may conclude that A is a maximal order. ■ 

By Corollary 12.71 we see that the fundamental parallelotope of the maximal order in Proposition 13.51 
has measure 10. Thus this code has 2.5 times the density of the Golden code. 

The algebra A3 has the drawback that the parameter 7 is quite large. This leads to an antenna power 
imbalance in both space and time domains. To some extent these problems can be alleviated by conjugating 
the matrix lattice by a suitable diagonal matrix (a trick used in at least [15]). One of the motifs underlying 
the perfect codes [10] is the requirement that the variable 7 should have a unit modulus. To meet this 
requirement we proceed to give a different construction for this algebra. 

Theorem 3.6: Let A be the square root of the complex number 2+i belonging to the first quadrant of the 
complex plane. The cyclic algebra QA+ = (Q(A)/Q(i), a, i), where the automorphism a is determined 
by <t(A) = —A, is a division algebra. The maximal orders of QA+ achieve the bound of Theorem 13.21 
Furthermore, the algebras QA+ and ^3 of Theorem 13.51 are isomorphic. 

Proof: The algebra QA+ is a central algebra F{u', A} over the field F = Q(z) defined by the 
relations A 2 = 2 + i, u' 2 = i, u'X = —An'. Comparing these relations with the relations in the proof of 

Theorem 13.51 we get an isomorphism of F-algebras / : QA-\ ► A3 by declaring f(u') = C /GV) = u 

and extending this in the natural way. The other claims follow immediately from this isomorphism and 
Theorem 13.51 ■ 

We refer to the algebra QA+ as the Golden+ algebra. This is partly motivated by the higher density 
and partly by the close relation between the algebra ^4 3 and the Golden algebra. After all, the algebra A3 
comes out when in the alternative description of the Golden algebra (cf. Example 12.11) the variable 7 = 5 
is replaced with its prime factor 2 + i. In Section |IV] we will provide an alternative proof for Theorem 13.61 
by explicitly producing a maximal order within QA+ and verifying that it has the prescribed discriminant. 
It is immediate from the discussion in the early parts of this section that in this case there is only one 
cyclic division algebra (up to isomorphism) with that discriminant. 

It turns out that all the algebras At in the Kiran-Rajan family of Example 12.51 have maximal orders 
achieving the discriminant bound. The following observation is the key to prove this. 

Lemma 3.7: Let F be either one of the fields Q(i) or Q(a>), and let Pi and P 2 be the two smallest 
ideals of its ring of integers R. Let V be a central division algebra over F, and let A be any i?-order in 
V. If the discriminant d(A) is divisible by no prime other than Pi and P 2 , then any maximal order T of 
V achieves the discriminant bound of Theorem 13.21 

Proof: We know that there exists a maximal order, say T containing A. The discriminant of T is 
then a factor of d(A), so Pi and P 2 are the only prime divisors of d(T ). From Theorem 13.11 we infer that 
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the only nontrivial Hasse invariants of V occur at Pi and P 2 . As the sum of the two Hasse invariants is 
an integer, they have the same denominator. This must then be equal to the index of V. The discriminant 
formula of Theorem 13.11 then shows that d(T ) equals the discriminant bound. Any other maximal order 
in V shares its discriminant with r . ■ 

Corollary 3.8: Let i > 2 be an integer. The maximal orders of the cyclic division algebra At = 
(Q(0) / QW> cr? 2 + i) from Example 12.51 achieve the discriminant bound. 

Proof: Proposition 12.51 and Lemma 12.91 indicate that the only prime factors of the discriminant of 
the natural order in At are 1 + i and 2 + i. The claim then follows from Lemma [3771 ■ 

At this point we remark that the natural orders of the algebras Ae of Example 12.51 are very far from 
being maximal. We will study this in greater detail in Section UV1 

Example 3.1: Let F = 3), so Of = Z[co}. In this case the two smallest prime ideals are generated 

by 2 and 1 — uj and they have norms 4 and 3, respectively. By Theorem 13.21 the minimal discriminant is 
4(1 — uj) 2 when n = 2. As the absolute value of 1 — uj is an application of the formula in Corollary 
12.81 shows that the lattice L of the code achieving this bound has m(L) = 27/4. In [27] we showed that 
a maximal order of the cyclic algebra (E/F, a(i) = —i, 7 = •>/— 3), where E = Q(i, V~ 3)> achieves this 
bound. 

We remark that one of the codes suggested in [15] is the natural order of the algebra of Example 
13.11 However, the authors there never mentioned the possibility of using a maximal order. Nor did they 
mention that their lattice actually is an order. 

IV. Finding maximal orders 

Consider again the family of cyclic division algebras At of index n = 2 l ~ 2 from Example 12.51 If is 
a maximal order of Ae, then according to Corollary 13.81 

d(A//Z[i]) = (1 + z) n(n_1) (2 + i) n{n ~ 1] . 
On the other hand, by Example 12.51 we know that 

d(At >nat /Z\i}) = (l + i)^ 2 M (2+i) n(r l -l)_ 

Hence, by Lemma 12.41 we may conclude that the natural order is of index 

[A, : A Anat ] = 2 «*- 5 > B+1 W 2 . 

In the cases t = 3, 4, 5 this index thus equals 2 3 , 2 26 , and 2 164 , respectively. In other words, using a 
maximal order as opposed to the natural order one can send 1.5, 6.5, or 20.5 more bits per channel use 
without compromising neither the transmission power nor the minimum determinant in the respective 
cases of 2, 4, or 8 antennas! Hence the problem of actually finding these maximal orders rather than 
simply knowing that they exist becomes quite relevant. In the following we shortly depict how maximal 
orders can be constructed in general. A more detailed version of the algorithm can be found in [18]. 

Let again F be an algebraic number field, A a finite dimensional central simple algebra over F, and A 
be a Z-order in A. Assume that A is given by relations (e.g. u 2 = 7), and that A is given by a Z-basis. For 
instance, we can always start with the natural order A (cf. Example 12.31) . We form a set S = {pi, ...,p r } 
consisting of the rational primes dividing d(A), i.e. A p is a maximal Z p -order if p S. 

The basic idea of the algorithm is to test for i = 1, ...,r whether A is maximal at p^. If the answer 
is yes, A is a maximal Z-order. If not, then at the first index % for which A p . is not maximal we can 
construct a Z-order r in A such that A p . C T Pi , and hence AcT (cf. Propostitions I2.10l - I2.15l ). This can 
basically be done in two steps. Let p E S. 

STEP 1 REPEAT UNTIL "YES": Compute I = 0- 1 (Rad(A p )) < A. Does the equality 0,(J) = A 
hold? 

"NO": Oi{X) D A 

A <- Oi(X) (Iteration step) 
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STEP 2 REPEAT UNTIL "NO": Compute the minimal ideals J u J 2 , ...,J h (h< dim Q A) of A/pA 
which contain Rad(A/pA). FOR i — 1, h compute Xj = (f)~ 1 (J'i). Does there exist an index i for which 
0,(Xi) > A? 

"YES": A «- C Z (X) (Iteration step) 
"NO": OUTPUT A is a maximal Z-order. 

Let p E S. First we test whether A p is an extremal (cf. Definition 12.71 ) Z p -order by checking if 
Oi(Rad(A p )) = Ap. If not, then we shall construct a Z-order T > A. If A p passes this test, then we can 
use the test of Proposition 12.151 If there exists an ideal J minimal among the ideals properly containing 
Rad(Ap) such that Oi(J) > A p , then we construct a Z-order T > A. Otherwise we correctly conclude 
that A is maximal at p and continue with the next p in the list S. In the end, the algorithm yields a 
Z-order A which is now maximal. The algorithm can be used similarly for constructing (9p-orders, but 
in the MAGMA software the implementations are for Z-orders only. 

For more details concerning the computation of the prime ideals in a ring, see [26]. A thorough 
explanation and an algorithm for computing the radical can be found in [28]. 

Let us next exemplify the above algorithm. 



A. 2 x 2 construction over Z[i] 

In the Golden division algebra (cf. Example 12. II or [3]), i.e. the cyclic algebra QA = (E/F, a, 7) gotten 
from the data E = Q(i, a/5), F = Q(z), 7 = i, n = 2, a(y/E) = —a/5, the natural order A is already 
maximal. The norm of the discriminant of A (with respect to Q) is 625, whereas the norm of the minimal 
discriminant is 100 [27]. We will now present a code constructed from a maximal order of the cyclic 
division algebra QA+ of Theorem 13 .61 The maximal order of QA+ also admits the minimal discriminant 
and is in that sense optimal. The algorithm now proceeds as follows. 

The natural order of the algebra QA+ is A = Z[i] © u'Z[i] © \Z[i] © u'XZ[i]. Hereafter, we will use 
a shorter notation A = (1, u 1 , X, u'X)zm for this. Let us consider A at the place P = 1 + i as it is the 
only factor of the discriminant for which we can enlarge A. The inverse image of the radical (|2.1 II) is 
J = 0" 1 (Rad(A/PA)) = ^((l + u',1 + X,l + u'X) Z2 ) = (1 + i,l + u',1 + X,l + u'X) m c A. A 
straightforward computation shows us (cf. Proposition 12.21) that the element 

l + u > + X + u'X (1 + M ')(1 + A) 

p ~ T+i " TTt e ° l{Jh 

which means that the answer to the question in Step 1 is "NO", and hence we set A' = (1, u', X, p)zu] 
and iterate. This time the inverse image of the radical is J' = _1 (Rad(A'/PA')) = + u',1 + 

X, 1 + p)z 2 ) = (l + i,l + u',l + X,l + p)z[i\ C A'. By taking the element 

r = ^ e O l{ J') 
1 + 1 

we can again enlarge the order A' to A" = (1, u 1 , r, p}zu] an d compute J" = 0~ 1 (Rad(A // /PA // )) = 
+ u', r, 1 + p)z 2 ) = (1 + h 1 + u',r, 1 + p)z[i] C A". We need one more iteration of Step 1. Now 
the element 

and the order A" is enlarged to A'" = (1, u, r, p)z[i\- From this iteration we finally get the answer to be 
"YES". 

In Step 2 there is nothing to do, as the only minimal ideal properly containing the radical is the radical 
itself. Hence we have constructed a maximal Z[i]-order of QA+ with a Z[z]-basis {1, u, r, p}. 

In order to give a concrete description of this order we describe it in terms of its Z[i] -basis. Let us 
again denote by A the first quadrant square root of 2 + i. The maximal order A consists of the matrices 
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aMi + bM 2 + cM 3 + dM 4 , where a, b, c, d are arbitrary Gaussian integers and M i; i = 1, 2, 3, 4 are the 
following matrices. 

M _ 1/(1-<)(1 + A) (l + i)(l-A)\ 
~~ 2 \(1 —«)(! + A) (1-0(1- A); • 



5. Enhancements to the Ivanyos-Rdnyai algorithm in some special cases 

The memory requirements of the above algorithm grow quite rapidly as a function of the dimension of 
the algebra. E.g. the MAGMA-implementation runs out of memory on a typical modern PC, when given 
the index 8 cyclic algebra A$ of Example 12.51 as an input. 

In this subsection we describe an algorithm that finds maximal orders for the algebras At. It is an 
adaptation of the Ivanyos-Ronyai algorithm that utilizes several facts special to this family of algebras. 
We list these simple facts in the following lemmas. We will denote Z[(i] by O for short. 

Lemma 4.1: The only prime ideal of O that lies above the prime 2 is the principal ideal P t generated 

by 1 - a 

Lemma 4.2: Let M be a finitely generated free C-module of rank k, and let and mi, . . . , m k be a basis. 
Let N be a submodule of M such that the index [M : N] is a power of two (in particular this index is 
finite). Then iV is also a free (9-module of rank k, and we can find a basis of iV of the form 

Hi = CLijlTlj , (Zjj E O. 

Proof: This is a straightforward modification of the proof of the corresponding result for modules 
over a PID. We briefly outline the argument, as we will need this later on. Let us start by choosing a basis 
mi,... , m k for M. We first consider the ideal I of those coefficients of m k that appear in expansions 
of elements of N. We have a natural surjective homomorphism from M/N onto O/I. Therefore the 
index of I in O is a power of two, so we may conclude that / is a power of the prime ideal Pi. By 
Lemma [47T1 1 is a principal ideal generated by a single element y k E O. We may thus choose an element 
n k = Vk m k + Si<fc %i m i from the submodule N. This will be the last element of a basis of N. We 
proceed by considering the submodule N' = N n J2i<k O m i °f vectors whose last coordinate vanishes. 
Then any element n E N can be written in the form n = z k nt + n' where n' E N' . The coefficients of 
m k ^\ that appear in N' then again form an ideal that by the Jordan-Holder theorem must be a power of 
Pi, and the argument can be repeated. In the end we get a free O -basis ni, n 2 , ■ ■ ■ , n k of N such that 

rii = ybjjmj, 
j 

where all the coefficients by E O, and by = whenever j > i. ■ 
Corollary 4.3: The maximal order Ag of At is a free (9-module of rank n = 2 e ~ 2 . 
Proof: We already know that A e contains A^ nat as a submodule of a finite index. Thus, there exists 
an integer M > with the property that MAe is a submodule of finite index in A^ nat . The formula for 
the discriminants tells us that we can further select the multiplier M to be a power of two. Clearly, it 
suffices to prove that MAg is a free module of the right rank. As the natural order, obviously, is a free 
(9-module of rank n, this is a consequence of Lemma 14.21 ■ 
Let then Y be any intermediate order, i.e. any order with the property A^ nat CTC A^. We will denote 
by T 2 the ring gotten by localizing Y at the prime 1 + i. This is naturally a subring of the corresponding 
localized version of the maximal order and consequently also of the completion of the maximal order Ai. 
This latter ring is a Z 2 [i]-order in the completion of the central simple Q2(0- ai g erjra Ai gotten from At 
by extending its scalars to the complete field Q 2 (0- Because the algebra At has a full local index 2 £ ~ 2 at 
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the prime 1 + i, At is actually a division algebra. By [20, Theorem 12.8] and the surrounding discussion 
therein we know that is a non-commutative discrete valuation ring, and that the (1 + z)-adic valuation 
of the reduced norm serves as a valuation. E.g. it yields a metric subject to the non-archimedean triangle 
inequality. So in the matrix representation the valuation of the determinant distinguishes the units from 
the non-units in the ring At. We immediately see that the same then holds in the ring T 2 — the units 
are precisely the elements whose reduced norm is a (1 + i)-adic unit. By the non-archimedean triangle 
inequality the non-units of T 2 then form its unique maximal ideal, which is then also the radical Rad(T 2 ). 

We summarize this line of reasoning in the following Lemma that is the key to our modifications to 
Step 1 in the main algorithm. 

Lemma 4.4: Let T be any intermediate order. The ideal X = T fl Rad(r 2 ) consists of exactly those 
matrices which determinants are divisible by 1 + i. 

The following lemma is a simple reformulation of the fact that P e is of index 2 in O. It will allow us 
to reduce the range of certain searches from O to the set {0, 1}. 

Lemma 4.5: Assume that p(x) = Yl^=oPi xt e Z[x]. Then 

p(Ce) = Po+Pi^ h Pk (mod P t ). 

Let us denote by s e the complex number 

^ = T ^ = 1 ^(i+o+o 2 + ---+cr 1 ). 

The fractional ideal generated by s e is then Pf 1 . 

Proposition 4.6: Let T be an intermediate order. Assume that it is a free O-module, and that gi, g 2 , . . . , g r , 
is its basis. Let / = </> _1 (Rad(r 2 )) (cf. Step 1). Then I is also a free (9-module of rank n that satisfies 
T C sil. We can find a basis for / that is of the form ri, r 2 , . . . , r n , where for all i either 

n = gi + ^2^ij9j, 

j<i 

such that all the coefficients G {0, 1}, or 

n = (i - Ct)g%- 

Proof: Any element of T has determinant (= its reduced norm) in Z[i]. The reduced norm of 1 — Q 
is an associate of 1 + i. Therefore (1 — Q)T C JCT. Thus the index of / in T is a power of two. Hence 
Lemma 14.21 implies that / is a free O-module of rank n. With the notation of Lemma 14.21 we also see 
that the coefficient y n is always either 1 or 1 — Q. In the former case Lemma [431 and the fact that 2 G Pi 
allow us to choose the coefficients as required. In the latter case we have no reason not to choose 
Ti — (1 — Ct)9i as this element is in / by Lemma I4~4l ■ 
Proposition 4.7: Let T, I, and the bases gi, . . . , g n and r 1; . . . ,r n be as in the previous proposition. 
Then the left order V = Ot(I) is a free O-module contained in SiT. It has a basis g[, . . . ,g' n , where for 
all i either 

g'i = Si(9i + ^2tij9j), 

j<i 

such that all the coefficients G {0, 1}, or 

9i = 9i- 

Proof: The inclusion (1 — Q)T C / immediately shows that T C Oi(I) C siT. Therefore the index 
of (1 — Q)V in T is a power of two. Again Lemma l4~2l shows that V is a free O-module. We also have 
the inclusion (1 — Q)V C T. An argument similar to the one in the proof of the previous proposition then 
shows that the algorithm in the proof of Lemma 14.21 yields a basis of the prescribed type. ■ 
When we use the natural order of the algebra At as a starting point, it is clear that p = 1 + i is the 
only interesting prime in Step 1 of the main algorithm. This step can now be completed simply by letting 
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T to be the natural order, and g 1 , . . . ,g n to be its (9-module basis. We next find a basis for Rad(T) by 
testing, whether any element of the type = g { + 2^2 j<i e ij9j has a determinant divisible by 1 + i (and 
if no such element is found then including r« = (1 — Q)9i into the basis instead). We then proceed to 
compute an (9-module basis for the left order V of this Rad(T). Again we simply check, whether any 
elements of the form g[ = s t (gi + 2~2j<i e ij9j) belong to V. Observe that it suffices to test a candidate of 
this form against the basis elements Vi only. If such an element is found, we record that V will be strictly 
larger than T. If no such element is found, we use g[ = gi instead. After we have done this for all i, we 
will know, whether V = T. If this is the case, we are done. Otherwise we replace T with V and repeat 
the process. 

We implemented this on the computer algebra system Mathematica, and on a typical modern PC it 
found a maximal order in the case I = 5 in less than half an hour. We believe that the memory savings 
due to the use of C-bases as opposed to Z-bases in the general purpose implementation in MAGMA 
account for this enhancement in the performance of the algorithm. This algorithm could naturally be 
ported into any CAS to handle these very specific cases. 

Example 4.1: Assume that we have the 4 antenna case £ = 4. Let us denote s = s<? for short. In this 
case the above algorithm yields an order with (left) (9-basis consisting of the elements ui, . . . ,w 4 : 

ui =1, 

u 2 =(s 2 + s 3 ) + s 3 u, 

u 3 =(s 4 + 2s 5 + 2s 6 + s 8 + s 10 ) + (s 5 + s 6 )u + s 10 u 2 , 

U4 =( S + s 4 + s 5 + s 8 + s 9 + s 10 + s ll + s 12 + s 13) + (g 9 + g ll + s 13 )u + (g 12 + § 13y + ^3 

We observe that the highest powers of s appearing in these basis elements are 0, 3, 10, and 13, respectively. 
This fits well together with our earlier calculation showing that the index of the natural order in a maximal 
one is 2 26 , as s _1 generates the prime ideal lying above 2, and + 3 + 10 + 13 = 26. 

It is a basic fact from the theory of the cyclotomic rings of integers that the conjugate of the element 
s is of the form a(s) = u a s, where u a is a unit of the ring Z[£]. Using this observation and the relation 
us = a(s)u we see that instead of the generator w 4 above we could use the product u 2 u 3 . After all, 
the (9-module spanned by these elements is an order, so we can utilize the fact that it is closed under 
multiplication. 

Example 4.2: In the 8 antenna case £ = 5 we get a free (9-module of rank 8 as a maximal order. The 
basis elements ui,...,ug are similar linear combinations of 1, u, w 2 , . . . , u 7 with coefficients of the form 
p(s), where p(x) E Z[x] and s = S£. In this case the polynomial coefficients of the various basis elements 
have maximal degrees 0, 3, 10, 13, 28, 31, 38, and 41. As expected, these degrees sum up to 164. Taking 
advantage of the fact that this module is also a ring we can describe the elements of the basis by 

ui =1, 

U2 =(S 2 + S 3 ) + S 3 U, 

u 3 =(s + s 2 + s 4 + 2s 5 + 2s 6 + s 8 + s 10 ) + (s 5 + s> + s w u 2 , 
ua =u 2 u 3 , 

s + 2s 2 + s 3 + 2s 4 + 5s 5 + 8s 6 + 8s 7 + 3s 8 + 5s 9 + 6s 10 + 5s 11 



u 5 



+ 7s 12 + 6s 13 + 7s 14 + 4s 15 + 5s 16 + 2s 18 + 2s 20 + g" + s< 



+ ( s 5 + 2s 6 + 4/ + s 8 + s 9 + s 10 + 2s ll + 2s 12 + 3s 13 + 3^4 + g 15 + 3^6^ 

+ (s 11 + 2s 14 + 2s 15 + s 16 + s 18 + s 20 )u 2 + (s 15 + s 16 )u 3 + s 2S u\ 

Uq =U 2 U 5 , 
U 7 =U 3 U 5 , 

u 8 =u 2 u 3 u 5 . 



18 



V. Analysis of the perfect algebras 

In this section we illustrate some computational techniques related to Hasse invariants and discriminants. 
We use the algebras underlying the perfect codes as test cases, because this may provide some additional 
insight into these algebras. This section places somewhat higher demands on the readers' background 
in algebra and algebraic number theory. It may be skipped if desired, as our code construction will not 
depend on the material in this section. 

Proposition 5.1: Let V\ = {E\/F, a%, 71) and T> 2 = (E 2 /F, a 2 , 72) be division algebras that have 
pairwise prime indices m\ and m 2 . Then T>\ ® V 2 is a division algebra with an index m\m 2 . Further, 

V X ®V 2 ~ (E^/F,^^,^ 2 ^ 1 ), 

where a x a 2 is an element of G&\{EiE 2 /F) ~ Gal(E 1 /F) x Gal(E 2 /F). 

Let Pi and P 2 be some pair of minimal prime ideals of the field F. If V 1 and V 2 have minimal 
discriminants that are only divisible by Pi and P 2 , then V\ ® V 2 has a minimal discriminant that is only 
divisible by Pi and P 2 . 

Proof: For the proof of the first two claims we refer the reader to [23, Theorem 20, p. 99]. The only 
nontrivial Hasse invariants of the division algebras T>i and T> 2 are those associated with primes P x and 
P 2 . The mappings in the fundamental exact sequence © are homomorphisms of groups. Together with 
the fact that extending scalars to a P-adic completion commutes with the formation of a tensor product 
shows that the Hasse invariants of T>\ V 2 are sums of those of T>\ and V 2 . Hence the discriminant of 
D\ ®V 2 is only divisible by the prime ideals P x and P 2 . By the proof of Theorem 13.21 it is then minimal. 

■ 

Suppose we have a finite cyclic extension E/F of algebraic number fields. Let P be a prime of F and 
B some prime of E that lies over P. We denote the completion E B by E P or E ■ F P . This notation is 
valid in Galois extensions, because the fields E B are isomorphic for all primes B that lie over P. 

A. 2 x 2 perfect code 

The first perfect algebra is the same as the Golden algebra QA = {E/F, a, 7), where the extension 
E/F = Q(i, v / 5)/Q(i) has discriminant (2 + i)(2 — i). The discriminant of the natural order is therefore 
(2 + z) 2 (2 — i) 2 . Because the discriminant of the algebra QA divides (2 + z) 2 (2 — i) 2 it can have at maximum 
two prime divisors (2 + i) and (2 — i). As a consequence the only Hasse invariants that can be nontrivial 

are h {2+i) and /i (2 -i). 

The algebra QA must have at least two nontrivial Hasse invariants and therefore hr 2+i \ and h( 2 -{\ are 
both nontrivial. Combining the equations LCM [m( 2 +j), ?7i(2-i)] = 2 and h/ 2+ i) + fo(2-i) = 1 we get that 
hr 2 +i) = h( 2 -{\ = 1/2. Theorem 13.11 states that the discriminant of QA is (2 + i) 2 (2 — i) 2 . Comparing this 
to the discriminant of the natural order we see that the natural order is maximal. 

B. 3 x 3 perfect code 

The underlying algebra of the 3x3 perfect code is P 3 = (E/F, a, oS), where again cu = (—1 + 3)/2, 
F — Q((v), E = Q(( 7 + Cf \ ^) and a : £7 + Cf 1 1 — > (? + Cf 2 - The algebra P 3 has a representation as 

L®u- L®u 2 ■ L 

where v? = uj. 

The discriminant of the extension E/F is (2 + 3) 2 (2 — v^— 3) 2 = Pf P 2 2 and the discriminant of the 
natural order has therefore only two prime factors. By Lemma 13771 the only nontrivial Hasse invariants of 
P 3 are h Pl and hp 2 . Because LCM [m Pl , m ft ] = 3. We get that mp 1 = mp 2 = 3. 

To calculate the Hasse invariant h Pl we pass to the completion Vp t = F Pl ®V 3 . From [20, Theorem 30.8] 
we get a cyclic representation Vp 1 = (E Pl /F Pl , a Pl , to), where E Pl /F Pl is a totally ramified extension 
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and a Pl is the natural extension of the automorphism a. Because the local index m Pl = 3, we know that 
Vp 1 is a division algebra. 

Next we try to find another cyclic representation for this algebra so that we can use the definition of 
Hasse invariant to calculate the value of hp,. 

It is readily verified that the field F Pl (u) = T Pl C Vp, is a cyclic and totally inert extension of F Pl . 
The Frobenius automorphism of the extension T Pl /F Pl is defined by the (Tp 1 / 'F Pl ,Pi)(u) = u 7 . The 
Noether-Skolem Theorem ([20, Theorem 7.21]) states that there is an element x E Vp 1 such that 

(Tp 1 /F Pl , Pi)(a) = x~ x ax Va G T Pl . (5) 

For an element x to fulfill © it is enough to satisfy the equation (Tp 1 /Fp 1 , -Pi) (it) = u 7 = xux^ 1 . By 
considering the equation ux = xu 7 = xto 2 u we see that x = Q 7 + Cf 1 + ^(C? + Cf 2 ) + ^((7 + Cf 4 ) £ L 
is a suitable element. 

We now prove that x 3 is an element of F Pl , and that vp 1 (x 3 ) = 1. The first statement follows from 
ua(x 3 ) = x 3 u = x 2 uuj 2 x = ux 3 . The second statement is obtained from the equation vp 1 (x 3 ) = 
v Pl (N E/F (x)) = v Pl (7(2 + (V=3» = L 

Proposition 16.41 now states that B\ = (T Pl /F Pl , (T Pl /F Pl , Pi), x 3 ) is a division algebra of index 3. By 
© we can consider B\ as a subset of the algebra V3. But B\ is a Fp x -central division algebra and hence a 
9 dimensional vector space over F Pl . From this we can conclude that (T Pl /F Pl , (T Pl /F Pl , Pi), a; 3 ) = V Pl . 

Lemma [6771 now implies that h Pl = 1/3. Because the sum of the Hasse invariants has to be an integer, 
the invariant hp 2 is 2/3. 

By considering the local indices we see that the discriminant of the maximal order is P 1 6 P 2 6 , that is, 
equal to the discriminant of the natural order. Thus, the natural order has to be maximal. 

C. 4 x 4 perfect code 

The division algebra under the 4 x 4 perfect code is V4 = (E/F, a, i), where Q(z) = F, Q(z, C15+C15 1 ) — 

£ and a : C15 + C15 1 — ► C15 + df- 

The extension E/Q(i) has discriminant d(E/Q(i)) = (2 + z) 3 (2 — z) 3 (3) 2 , and the only Hasse invariants 
that can be nontrivial are h^,h^ 2+ i) and hi^-i). We use similar methods to those in the case of P 3 to 
get that h( 2 +i) = 3/4 and h(z-i) = 1/4- The sum h( 2 -i) + = 1 an d therefore h(s) must be trivial. 

Further, the local indices reveal that the discriminant of the algebra is (2 + i) 12 (2 — i) 12 . The discriminant 
of the natural order on the other hand is (2 + i) 12 (2 — i) 12 (3) 8 . Lemma 12741 tells us that the index of the 
natural order in the maximal order is 81. 

D. 6 x 6 perfect code 

In the 6x6 perfect code construction the center is F = Q(uj) and the maximal subfield E = K(9), 
where 9 = ( 2 s + C^- 

In [10] where the perfect codes were introduced, the authors gave the mapping a\ by the equation 
°i : C28 + C28 1 1 — > C28 + C28 2 - Unfortunately, this mapping is not an P-automorphism of the field E. We 
replace <j\ with the automorphism a defined by the equation a : C28 + C28 1 1 — * Cls + C28 5 - The relative 
discriminant of the extension E/F is (2) 6 (2 + v /= 3) 5 (2 - v^) 5 = (2) 6 (7) 5 . We denote the resulting 
algebra by V§. 

Thus the Hasse invariants of Vq that can be nontrivial are h^ 2+y ^3), h/ 2 -^z^, and h^y 

Now we are going to present V% as a product of two smaller division algebras. We first calculate the 
Hasse invariants of these smaller algebras and then from these derive the Hasse invariants of V§. 

Let us first consider the algebra B 2 = (Q(v7, u)/Q(u), a 2 , —uS). The algebra B 2 is a division algebra 
with Hasse invariants /i( 2 - v /^3) — ^(2+v^3) = 1/2- The proof is postponed until the end of Section [VTl 

The algebra V3 = (E/F,a 3 ,u) was previously shown to be a division algebra with Hasse invariants 
^(2- v /3 3) = 2/3 and /^( 2 + v ^3) = 1/3- We now consider the algebra B 3 = (E / P, cr 3 , to 2 ) . By [20, Theorem 
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30.4] we have V 3 <g> £> 3 ~ (E/F,a 3 , 1) ~ M 3 (F). This shows that V 3 ® i3 3 has trivial Hasse invariants 
and therefore the Hasse invariants of B 3 are fy 2 _^z3) = 1/3 and hr 2+y /=s\ = 2/3. 
If we now consider the algebra B 3 B 2 = 

(Q(V7, w) ■ Q(C 7 + C 7 ~\ w)/Q(w), a 2 a 3 , (-^) 3 ■ (to 2 ) 2 ) 

it is seen that the corresponding Hasse invariants are /i( 2 _ v /r 3 ) = 1/3 + 1/2 = 5/6 and h 2+ ^^ = 
1/2 + 2/3 = 1/6 (mod 1). 

By considering the equation cr 3 (C 7 + Cf 1 ) = C? + Cf 2 = Cl + Cf 5 we notice that (720-3 = (76. Combining 
this and the equation (— cu) 3 • = — uo we get that £> 3 (g> £> 2 ~ V®. 

The algebra P 6 has only two nontrivial Hasse invariants that are = 5/6 and fy 2 -V=3) = V^- 

Whence, the discriminant of the maximal order is (2 — -\/— 3) 30 (2 + V— 3) 30 = (7) 30 . The discriminant of 
the natural order on the other hand is (2) 36 ( 7) 30 . In this case Lemma [2~4l tells us that the perfect lattice 
is of relatively high index 2 18 in its counterpart within the maximal order of same minimum determinant. 



VI. Constructing division algebras with a minimal discriminant 

We have divided this section into two parts. In the first part we are concentrating on algebras that have 
a cyclic representation with a unit non-norm element 7. 

In the second section we relax the restriction on the size of 7 and we give a general construction for 
Q(z) and Q(\/— 3) -central division algebras with a minimal discriminant. 

One should note that none of the natural orders of the algebras we shall construct has a minimal 
discriminant. This, unfortunately, is not just a coincidence. In the following we prove that there are no 
natural orders reaching the bound of Theorem 13.21 

In the next lemma we use some basic results from the theory of discriminants and differents. For these 
results and the notion of different we refer the reader to [29, Chapter 3.12]. 

Lemma 6.1: Suppose we have a Galois extension E/F of degree n and that there are g prime ideals 
E>i of E lying over the prime P of F. If the prime P is wildly ramified in the extension E/F, then 

v P (d(E/F)) > n. 

Proof: Suppose that De/f is the different of the extension E/F. Then it is an easy exercise in 
Galois theory to show that vbXDe/f) — v Bj(D E / F ) for every i and j. Because P was supposed to be 
wildly ramified 

s = v Bi (D e /f) > e, (6) 

where e is the ramification index of Bi/P. 

The theory of normal extension states that efg = n, where / is the inertial degree of Bi/P. Taking 
into account this and © we can conclude that 

v P {d{E/F)) = v P (N E/F (D E/F )) = sgf > egf = n. 

■ 

Proposition 6.2: Suppose we have a division algebra V = (E/Q(i), a, 7), where E/Q(i) = n and 7 
is an algebraic integer. If A is the natural order of the division algebra V, then 

|d(A/O Q(0 )| > |(2 + z)"("- 1 ))(l+i)^- 1 )|. 
Proof: The natural order A is a subset of some maximal order A max and therefore |<i(A/C , Q(i))| > 
|(2 + i)"("- 1 )(l + i)"("- 1 )|. Let us then assume that \d(A/0 Q{i) )\ = |(2 + i)"( n ^ 1 )(l + i) n{ - n -% 

According to Lemma [2791 the only primes that could be ramified in the extension E/Q(i) are (1 + i), 
(2 + i), and (2 — i). Lemma [67T1 assures that none of these primes could be wildly ramified. 

One of the main results of the global class field theory [24, p. 124] states that there exists a ray class 
field C(i + j)( 2 +j)(2-i) that contains all the cyclic extensions of Q(z) where (1 + i), (2 + i), or (2 — i) is 
tamely ramified. 
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We can now calculate the degree of the extension C( 2 +i)(i+i)(2-i)/Q(i)- By [24, Theorem 1.5] we have 
[C (2+ i)(i + i)(2-i) : Q(*)] = 2 > which implies that E = C(2+t)(i+i)(2-i) and n = 2. 

The ray class fields C/2+i)(i+i) and C(2-i)n+i) that admit tame ramification at (2 + i) and (1 + i) or, at 
(2 + z) and (1 — i), respectively, are both trivial extensions of Q(i). Hence, both (2 + i) and (2 — i) are 
ramified in E and divide the discriminant of the extension E/Q(i), The discriminant of the natural order 
A now has to be divisible by at least (2 + i) 2 (2 — i) 2 . This gives us a contradiction. ■ 

Proposition 6.3: Suppose we have a division algebra V = (E/Q(y/—3), a, 7), where E/Q(y/—3) = n 
and 7 is an algebraic integer. If A is the natural order of the division algebra V, then 

|d(A/0 Q( v=S))l > \(V=S) n{n - 1)) (2r^\. 
Proof: The proof is similar to that of the previous proposition. ■ 
These considerations reveal that reaching the optimal density of a code-lattice requires considering 
maximal orders instead of natural ones. 

We give one simple lemma for later use, it is a slight generalization to [11, Theorem 1]. We denote 
the multiplicative ideal group of the field F by (If)*- 

Lemma 6.4: Let E be a Galois extension of a number field F and let P be a prime ideal of Of that 
lies under the prime B of the ring F - If the inertial degree of P in the extension E/F is / and 7 is 
such an element of F that (vp(j), f) = 1, then Y Ne/f(E) for any i = 1, 2, ...,/— 1. 

Proof: The ideal norm of B is N e /f(B) = where / is the inertial degree of P in the 
extension E/F. It is clear that the group N e /f((If)*) is generated by the norms of prime ideals and that 
{N E/F (a)0 F \aeE*} C N E/F (I F ). Therefore f\v P (N E/F (a)0 F ) for all a E E. ■ 

A. Algebras with a unit 7 

TABLE I 

Q(i)-CENTRAL DIVISION ALGEBRAS WITH A UNIT 7 



n 


7 


/« 


2 


i 


a; 2 - 


-(2 + i) 


4 


i 




"(2 + 



ij Center Q(i): In Table U we give a cyclic representations for algebras of degree 2 and 4. Proposition 
12. II implies that 4 is the biggest degree that we can hope to have a cyclic division algebra with a unit 7. 
There does not exist such an algebra of degree 3. The reason for this is that in every cyclic extension 
E/Q(i) of degree three, all the units of Q(i) are third powers and therefore are in the image of the norm 

N E /Q(i). 

In the following we use the generic notation Q(i) = F and E = F(a n ), where a n is a zero of the 
polynomial f n . 

Algebra V 2 : The algebra V 2 was previously shown to be a division algebra with a minimal discriminant. 

Algebra V A : When considering V A we first have to check whether it really is a division algebra. We 
note that (2 + i) is a totally ramified prime in E/F. This results in the local extension E(2+i)/F^+i) being 
a totally and tamely ramified cyclic extension of degree 4. We note that #(O f . 2+ ../(2 + i)0 F(2+1 ) = 
4(0 F /(2 + i)) =5. 

Proposition 12.11 states that £> 4 is a division algebra if i satisfies the norm condition, i.e. neither of the 
elements {i, — 1} is a norm. 

Hasse Norm Theorem [20, Theorem 32.8] states that it is enough to show that the elements {i, — 1} 
are not norms in the extension E( 2 +i)/F( 2 +i). Elementary local theory [30, Proposition 7.19] states that if 
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we have any complete residue system {0, 1, a, b, c} of the group /(2 + and an arbitrary 

unit e e P(2+i) then 

E* +i = {1, a, b, c} x (1 + (2 + i)^) x (e(2 + i)). (7) 

The prime (2 + i) is tamely ramified in Er 2+i \ /P(2-b) an d therefore the local conductor is (2 + i) ([24, p. 
12]). The definition of the conductor now implies that (1 + (2 +i)0^ 2+ . C iVg 2+ . /p 2+ . {E(2+i))- Because 
the prime (2 + i) is totally ramified, we have ei(2 + z) C ,^ (-F} 2+ i)) for some unit ei e P(2+i)- 

The previous results now imply that (1 + (2 + i)^,^.,) x (ei(2 + i)) C N$ /P (E 2+i ). 

One of the main theorems of local class field theory states that (F( 2 +i))* /(Ng^ + ■ ) /P 2+ ■ = 
Gal(Er2+i)/F(2+i))- By considering © we see that the elements {a,b,c} are not norms. Because the 
elements {0, i, — 1, —i, 1} form a complete residue system of the group Op {2+ /(2 + z)C?^ + we find that 
neither of the elements {i, —1} is a norm. 

The discriminant of the extension E/F has only two prime divisors (2+i) and (1+i) and therefore also 
the discriminant of the natural order of P4 has only two prime divisors. This implies that the discriminant 
of the algebra is minimal. 

TABLE II 

Q(c<j) -CENTRAL DIVISION ALGEBRAS WITH A UNIT 7 



n 


7 


fn 


2 


—uj 


x 2 


+ V=3 




3 


LU 


x 3 


- 2 




6 


-LU 2 


x 6 


- ?>yf^ix A + 4x 3 - 9x 2 + 12V 


-3z + 3 v /Z 3 + 4 



2) Center Q(\/— 3): In Table [XT] we give cyclic representations for algebras of degrees 2, 3, and 6. The 
theorem of Albert shows that 6 is the biggest degree we could hope to have a division algebra with a unit 
7. We cannot have a division algebras of degrees 4 and 5 as tensoring these with a division algebra Q 3 
(below) would respectively give us division algebras of degrees 12 and 15 with a unit 7. 

We use the same generic notation as in the case of Q(i)-central algebras. 

Algebra Q 2 : We use here the same methods that were used with the algebra V^. We remark that 
(V - 3) = P is tamely ramified in the extension E/F. If we pass to the completion Ep/Fp we get that 
the local conductor is P and that {—to, 1,0} is a complete set of representatives of the group Op P /p- As 
a result it is seen that —to is not a norm in the extension Ep/Fp and therefore it is not a norm in the 
extension E/F either. From this it follows that Q 2 is a division algebra. 

By now it is obvious that the discriminant of the natural order of the algebra Q 2 has only two divisors 
(yf— 3), and (2) and hence the maximal order admits a minimal discriminant. 

Algebra The proof of this case is similar to that of Q 2 except that the tamely ramified prime P is 
2 and that the suitable set of representatives is {1,lj,lj 2 }. 

Algebra Q 6 : The algebra Q 6 we got as a tensor product from the algebras Q 2 and Q 3 . 

The postponed proof. When we were discussing the 6x6 perfect code we postponed the analysis of 
the algebra B 2 = (E/F, a 2 , —to), where E/F = Q(\/7,j)/Q(u). Now we have enough methods to attack 
this problem. We use similar strategy as in the case of the algebra V4. 

The prime (2 + V - 3) = Pi is tamely ramified in the extension E/F. By passing to the Pi-adic 
completion E Pl /F Pl we find that the local conductor is P 1 . The image of the norm N^ p ,p p includes 

((1 + Pi)) x (e(2 + V^3)), where e is a unit of F Pl . 
The set {0, l,u, —uj,uj 2 , —uj 2 } is a complete residue system of the group Op p j 'P\Op p ^ and whence 

(P Pi r = (-j)x(l + Pi)x(e(2 + v / =3)). 
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On the other hand #((F Pl )* / N Ep /p (EpJ) = 2 and therefore —j cannot be a norm. From this it follows 
that the local algebra (B 2 ) p 1 is a division algebra of index two. 
There is no other choice for the Hasse invariant h Pl than 1/2. 

Replacing the prime P 1 with P 2 = (2 — \/— 3) in previous considerations we see that h P , 2 = 1/2. 

The extension E/F has only three ramified primes (2— (2+v^— 3), and (2). Thus, the discriminant 
of the algebra i3 2 can have three prime divisors at maximum. The potential nontrivial Hasse invariants of 
£> 2 are now h Pl , hp 2 , and ht 2 y The sum of h Pl and hp 2 is 1 and therefore ht 2 \ must be trivial. 

B. General construction 

In their recent paper [14] Elia et al. gave an explicit construction for division algebras of an arbitrary 
degree with centers Q(i) and Q(\/— 3). In their general constructions they used non-unit, but relatively 
small 7's. As they were not interested in maximal orders nor the discriminants of the corresponding 
division algebras their algebras (with few exceptions) did not happen to have minimal discriminants. 

We are now going to give a general construction for division algebras of arbitrary degree and with 
minimal discriminants. Due to Proposition 15.11 we can concentrate on algebras of prime power index. As 
a drawback our constructions will be dependent on the existence of certain prime numbers. We discuss 
this existence problem in Section IVI-CI which is purely number theoretic. 

We first consider two easy prime powers and then move forward to more complicated ones. 

For ease of notation in this subsection we will denote by Z m the residue class ring modulo m, i.e. 
Z m = Z/mZ. Thus e.g. is logically the group of units of that ring. 

Lemma 6.5: Suppose that E is a cyclic extension of F and that aOp = Pi and P 2 are a pair of smallest 
primes in F. Assume that P\ is totally inert and P 2 is the only ramified prime in the extension E/F. 
Then 

A — (E/F, a, a), 

where (a) = G&\(E/F), is a division algebra that has a minimal discriminant. 

Proof: Lemma |6~4l combined with Proposition 12. 1 1 gives that A is a division algebra. The minimality 
of the discriminant follows from Lemma 13.71 ■ 

Example 6.1: Lemma [631 is nothing but a simple generalization of Corollary 13.81 where we gave a 
construction for a family of Q(z) -central division algebras of degree 2 k with a minimal discriminant. 

Example 6.2: The field Q(C3 fc +0 has a unique subfield Z with [Z : Q] = 3 fc . The extension Q(y/—3)Z/ 
Q(V— 3) has degree 3 fe and the prime (2) is totally inert in this extension. The extension also has a very 
limited ramification, the prime ( v / — 3) is the only ramified one. 

Primes (y/— 3) and (2) are a pair of minimal primes in the field Q(\/— 3). Lemma 1631 states now that 
the cyclic algebra A = (Q(v^— 3)Z/Q(y/— 3), a, 2) is a division algebra with a minimal discriminant. 

In Example 16.21 we found a suitable extension £'/Q(v^— 3) that only had one ramified prime (V— 3). 
However we can prove that for an arbitrary degree there usually does not exist a cyclic extension that has 
ramification over (V— 3) or (2) only. This assures that in general we cannot use such simple methods. 
Next we will provide a construction method that takes care of most of the prime power degrees. First we 
need some preliminary results. 

We now present a global Frobenius automorphism. Suppose we have a finite Galois extension E/F 
and that B is such a prime ideal of Op that B n Op = P is unramified in the extension E/F. There 
exists a unique element (B, E/F) of the group Gal(E/F) that is associated to the prime B. We call this 
element the Frobenius automorphism of B. 

If the extension E/F is abelian, all the primes Bi that lie over P have the same Frobenius automorphism 
and we can denote (B, E/F) by (P,E/F). 

For the definition and properties of the Frobenius automorphism we refer the reader to [31, p. 379]. 

We consider a tower of fields F\ C F 2 C E of finite extensions. 
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Proposition 6.6: If F± Q F 2 Q E, E / Fi and P 2 /Pi are normal and P is such a prime ideal of E that 
P fl Fi = P is unramified in EjF-y, then 

{B,E/F X )\ F2 = (BnF 2 ,F 2 /F 1 ). 

The prime P is totally inert in the extension E/F\ if and only if (P, E/Fx) generates the group Gal(P/Pi). 
Proof: [31, Theorem 7.10, p. 380]. ■ 
The next lemma is a rather direct consequence of the definition of Hasse invariant. 
Lemma 6.7: Let 

A=(E/F,a, 1 ) 

be a division algebra where (a) = G(E/F), 7 £ F*, [E : F] — n and suppose that P is a prime ideal of 
F that is totally inert in the extension E/F. If k is the smallest possible positive integer so that a k is the 
Frobenius automorphism of P then the Hasse invariant of P 

tip = . 

Proof: [20, p. 281]. " ■ 

Let us next consider a tower of fields Pi C P 2 C P of finite extensions and the proofs of the next two 
simple lemmas will be omitted. 

Lemma 6.8: Let B be a prime ideal of E, P 2 = Of 2 fl B and P\ — Op x D P. 

1. Let f(B/Pi), f(B/P 2 ), and f{P 2 /P\) be the respective inertia degrees of P over Pi, P over P 2 , 
and P 9 over P 1 . Then 

/(P/PO = f(B/P 2 )f(P 2 /P 1 ). 

2. Let e(P/Pi), e(B/P 2 ), and e(P 2 /Pi) be the respective ramification indices of P over Pi, P over 
P 2 , and P 2 over P x . Then 

e(P/Pa) = e(P/P 2 )e(P 2 /P 1 ). 
Lemma 6.9: Let E/F be a Galois extension, P a prime ideal of P and P = F f] B. Then 

e(P/P) I [P : P] 

and 

/(P/P) I [P:P]. 

Lemma 6.10: Let p be a prime and n such an integer that n\(p — 1). The field Q(C P ) has a unique 
subfield Z with [Z : Q] = n. 

There exists a group isomorphism from Z*/(Z*) ra to Gal(Z/Q) that takes any prime pt 7^ p to the 
corresponding Frobenius automorphism (pi,Z/Q) in Gal(Z/Q). 

The prime pi 7^ p is totally inert in the extension Z/Q if and only if p* is not an nth power (mod p) 
for £ = 1, . . . , n — 1. 

Proof: It is well known that there exists a unique isomorphism ip from Z* to Gal(Q(£ p )/Q) which 
takes prime p\ 7^ p to (pi, Q(Cp)/Q)- We denote the fixed field of the group ^(Z*) n by Z. It is now 
clear that Z is unique and [Z : Q] = n. If we first map the elements of Z* with ip to Gal(Q(C P )/Q) and 
then restrict the resulting automorphisms to the field Z, we obtain an isomorphism <p from Z*/ (Z*) n to 
Gal(Z/Q). Proposition 16.61 states that has the claimed properties. 

The last claim follows from the properties of cf> combined with the last statement of Proposition |6.6i ■ 
Proposition 6.11: Suppose that P = Q(y / c) is a quadratic field, q ^ 2 is a given prime and n a given 
integer. We suppose that P 1 and P 2 are the smallest primes ideals in P and pi and p 2 are the prime 
numbers that lie under Pi and P 2 . 

Let p be such a prime that q n \(p — 1), (p, c) = 1 and that pi and p 2 are totally inert in the extension 
Z/Q, where Z is the unique subfield of Q(( p ) of degree q n . We also suppose that p is inert in the 
extension P/Q. 
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The extension FZ/F is a cyclic Galois extension of degree q n where the prime ideals Pi and P 2 are 
totally inert and P = pO F is the only ramified prime ideal in the extension FZ/F. 

Proof: Let B be a prime ideal of FZ, P z = O z n B, P F = O f n B and b = Q n B. We denote 
the corresponding ramification indices by e(B/P z ), e[P z /P F ) and e(P F /b). According to Lemma [6781 

e(P/6) = e(B/P z )e(P z /b) = e(B/P F )e(P F /b). 

Lemma [6791 for its part states that e(B/P z ),e(P F /b) | 2 and e{P z /b),e{B/P F ) \ q n . This together with 
the previous equation shows that the prime P F C O f is ramified in the extension FZ/F if and only if 
the prime 6 is ramified in the extension Z/Q. 

The prime p is the only ramified prime in Z/Q and because p is inert in the extension F/Q we see 
that P is the only ramified ideal in the extension ZF/F. 

If we choose B so that P F = Pi or P F = P 2 , then 

f(B/b) = f(B/P z )f(P z /b) = f(B/P F )f(P F /b) = q n -c } 

where c = 1 or c = 2. This combined with Lemma |679l implies that f(B/P F ) = q n . ■ 
In the following propositions we use the notation from Proposition 16.1 II 

Proposition 6.12: There exists such a group isomorphism between Gal(FZ/F) and Gal(Z/Q) that 
every Frobenius automorphism of B C O fz maps to the Frobenius automorphism of B PI Z = B z . 

Proof: It is a well-known fact that there exists a well defined surjective homomorphism from 
Gal(FZ/Q) to Gal(Z/Q) for which o \ — > o\ z . The kernel of this map consists of those elements of 
Gal(PZ/ Q) that act trivially on the field Z. On the other hand, if we restrict the domain of the map to 
those elements that act trivially on F this map is an injection because the only element of Gal(FZ/Q) that 
acts trivially on both fields F and Z is the identity map. As we know that \Gal(FZ/F)\ = |Gal(Z/Q)| 
the described map must be an isomorphism. Now the statement about Frobenius maps follows from 
Proposition 16.61 ■ 

Proposition 6.13: Let 

P2P1 = 1 (8) 

in the group Z* p /{Z* p )i n , Pi = a x Q F , and P 2 = a 2 F . Then 

A = (FZ/F, a, oia a ) 

with (a) = Gal(FZ/F) is a division algebra that has a minimal discriminant. 

Proof: The prime P 1 is totally inert in the extension FZ/F. Thus, Lemma [6741 states that A is a 
division algebra. 

From the cyclic presentation of the algebra A we instantly see that A has only three Hasse invariants 
that can be nontrivial: h Pl , hp 2 , and hp. In what follows we are going to show that the invariant hp must 
be trivial. 

We first choose o to be the Frobenius automorphism of Pi. Lemma [6771 now shows that the Hasse 
invariant of Pi is 

qn 

Because the group Z*/ (Z*) 9 ™ is cyclic we get from ® that p 2 = pf~ l in Z*/ (Z*) 9 ". This implies that 
(P 2 ,FZ/F) = a n - 1 . Lemma E3 then states that 

q n - 1 



q n 



h 



p 2 . 



The sum of the Hasse invariants of A must be zero (mod 1), whence 



h Pl + hp 2 + hp G Z. 
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But, we already saw that h Pl + hp 2 £ Z, which implies that hp £ Z. The discriminant of the algebra *4 
has now only two divisors Pi and P 2 . 

In the beginning of our proof we make the assumption that a is the Frobenius of the prime P 1 . 
However, the choice of the generator of the group Gal(FZ/F) in a cyclic representation does not change 
the discriminant of the corresponding algebra. ■ 

Example 6.3: Suppose that the center F = Q(i). The primes (1 + i) and (2 + i) are a pair of smallest 
prime ideals in this field. We want to produce a division algebra of index 10 that has a minimal discriminant. 
It is not difficult to check that 2 t and 5* are not 5th powers (mod 11) for t = 1, ... ,4, and that 11 is 
inert in the extension F/Q. Lemma 16. Ill states that Q(Cu) has a subfield Z, [Z : Q] = 5, and that 2 and 
5 are totally inert in the extension Z/Q. 

Proposition 16. 1 1 1 states that the primes (1 + i) and (2 + i) are totally inert in the extension FZ/F and 
the prime ideal HOp is the only ramified ideal in the extension FZ/F. 

We easily see that 2 ■ 5 = 1 in Z* 1 /(Z* 1 ) 5 . Therefore, 

(FZ/F, o-i, (1+i) (2 + i)) 

is a division algebra with a minimal discriminant. 

We previously saw that A = (Q(C2 4 )/F, a 2 , 2 + i) is a division algebra of index 2 and has a minimal 
discriminant. Finally, from Proposition 15.11 

(Q(C 2 4)Z/F, ( x 1 a 2 ,(l + z) 2 (2 + z) 7 ) 

is seen to be a division algebra of degree 10 with a minimal discriminant. 



C. Existence of suitable primes 

Propositions 16.1 II and 16. 131 have turned our construction project into a hunt of suitable prime numbers. 
The problem is that we do not know if there are "enough" suitable prime numbers. The answer is that in 
most cases there are. This will be proved in Theorem 16.171 but first we need some preliminary results. 

For the definition and the basic properties of Kummer extensions we refer the reader to [29, p. 197]. 

Proposition 6.14: Let E/F be a Kummer extension with E = F(a), a n = a £ Op, and let P be a 
prime ideal of F that is not a divisor of a ■ n. Furthermore, let t be the largest divisor of n such that the 
congruence 

x l = a (mod p) 

has a solution in Op. Then P decomposes in E into a product of t prime ideals of degree n/t over P. 
Proof: [29, Theorem 6.8.4, p. 197]. ■ 
Lemma 6.15: Suppose that q and p are prime numbers and that q t \(p — 1) for some integer t. If c is 
an integer and the equation 

c = x q (mod p) (9) 
is not solvable, then neither is any of the equations 

c k = x qt (modp), (10) 

where k = 1, . . . , q* — 1. 

Proof: Let a be a generator of the cyclic group Z* Then we can write that c = a n (mod p) for 
some integer n. 

Let us assume that © has no solution. This implies that q is not a factor of n. Assume then that for 
some k there is a solution d for (flOl) . If we write d = a s , then (flOl) gives that kn — sq l = v (p — 1), where 
v is some integer. As q t \(p — 1) this would mean that q t \kn. That gives us a contradiction. ■ 

In the following we use the phrase "the prime P has inertia in the extension E/F". By that we mean 
that at least one prime ideal B of E that lies over the P has inertial degree f(P\B) > 1. 
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Lemma 6.16: Suppose that Fx and F 2 are Galois extensions of a field F and F\ fl F 2 = F. The prime 
P of Op has inertia in the extension FiF 2 if and only if it has inertia in the extension Fx or F 2 . The 
prime P is ramified in the extension FxF 2 if and only if it is ramified in Fx or in F 2 . 

Proof: For the proof the reader is referred to [32, p. 263]. ■ 

The proof of the following theorem is a slightly modified version of the proof of [33, Theorem 1]. 
We do not suppose here that the center is totally complex nor that the ring Op is a PID. However, we 
suppose that p x 7^ Vi- 

Theorem 6.17: Assume that F = Q(y / c) is a quadratic field, Pi and P 2 are the smallest primes in F, 
q 2 is a given prime, and n a given integer. Let us also suppose that p x and p 2 axe prime numbers that 
lie under Px and P 2 . 

If q \ c, then there exists infinitely many prime numbers p so that p is inert in F, Q(Cp) has a unique 
subfield Z, [Z : Q] = q n , where px and p 2 are totally inert, and pxp 2 = 1 in Z*/(Z*) 9 ™. 

iVoo/": Let us denote q n = s, Q(C)((PiP2) 1/s ) = K, K((px) l/q ) = K x and suppose that q ^ px. By 
considering the prime ideal factorization of pip 2 in Q(C«) we mav conclude that (pxP2) d cannot be an sth 
power for any d — 1, . . . , s — 1. Therefore [if : Q(Cs)] = 

As we have supposed that q \ c there has to be at least one prime p 3 that has a ramification index 2 in 
the extension F/Q, but is not ramified in the extension Q(Cs)/Q- Earlier, we saw that [K : Q(Cs)] — s - 
Because p 3 is not ramified in F/Q and 2 does not divide [K : Q(Cs)]* none of the prime ideals P 3 in Ok 
that lies over p 3 has 2 as a divisor of the ramification index e(P 3 \p 3 ). This implies that F K . 

By [33, Lemma 2] we know that [Kx : K) — q. Because q ^ 2 and F % K the extension KxF/K is 
cyclic and [ifiF : if] = 2g. 

Chebotarev's density theorem [31, Lemma 7.14, p. 392] states that K has infinitely many prime ideals 
that have absolute degree one and are totally inert in the extension K{^/p{)F/ K. We choose one, P, that 
not only has an absolute degree one but that is also unramified in the extension K/Q. 

We denote the prime of Q that lies under P by p. The field Q(Cg n ) is a subfield of K and therefore p 
splits completely in the extension Q(Cg»)/Q- The theory of cyclotomic fields [29, p. 195] now gives that 

p = 1 (mod q n ). 

Next we are going to show that p\ is not an sth power (mod p) for t — 0, . . . , s — 1. We assume the 
contrary. Suppose that px = a q (mod p) for some integer a. Now px = a q (mod P). This last equation 
however cannot be true because P is totally inert in the Kummer extension Kx/K. Lemma [6.151 now 
states that equation px = x % (mod p) does not have a solution for any t = 1, . . . , q n — 1. 

Lemma 16.101 states that Q(C P ) has a unique subfield Z with [Z : Q] = g n , and that px is totally inert 
in the extension Z/Q. 

The prime P has absolute degree one in K and therefore (pip 2 ) 1 ^™ = c (mod P), where c is some 
integer. This implies that 

Pxp 2 = c q (mod p). 

If we use the notation of Lemma 16 .10[ the map takes px to the generator g of the group Gal(Z/Q) and 
Pi ■ p 2 to identity. The map is a homomorphism and therefore 4>{p 2 ) = g~ x , which again is a generator 
of the group Gal(Z/Q). Lemma [6. 101 now shows that p 2 is totally inert in the extension Z/Q. 

To complete the proof we have to show that the prime p is inert in the extension F/Q. The prime P 
must be inert in the extension FK/K and therefore the prime p has at least some inertia in the extension 
FK/Q. Because p is totally split in the extension K/Q it does not have any inertia in this extension and 
therefore Lemma 16.161 states that p must be inert in the extension F/Q. ■ 

Theorem 16. 171 states that for the center Q(i) the only problematic prime power indices are of the form 
2 k . Luckily, the construction of Corollary 13.81 covers these cases. As a consequence, we can construct 
a division algebra with a minimal discriminant for an arbitrary index. In Table [III] we give explicit 
representations for division algebras with a prime power index (< 20) and a minimal discriminant. 



28 



For each index q n we have searched the prime p of the Theorem 16.171 along the lines of Example 16.31 
After the prime p is found the actual minimal polynomial of the extension FZ/FQ can be easily found 
by considering the subfields of the extension Q(C P )/Q. Both tasks were done with the aid of computer 
algebra system PARI [34]. 

If the center is Q(\/— 3), the problematic prime powers are 2 n and 3 n . Algebras of degree 3" we get 
from Example 16.21 but degrees 2 n are more problematic. Still for indices 2 and 4 we can find suitable 
primes even when Theorem 16. 171 is not promising anything. As a conclusion we can construct a division 
algebra with a minimal discriminant if the index is not divisible by 8. 

In Table [IV] we give explicit representations for our algebras. 

Example 6.4: From Table [III] we get that 

^3 = (Q(i)(a 3 )/Q(i),«T 3 ,(l+i)(2 + i)) 

and 

A 2 = (Q(i)(a 2 )/Q(i),a 2 ,(2 + i)) 

are division algebras with minimal discriminants. According to Proposition 15.11 algebra A 2 <8> A3 = 

(Q(i)(a 6 )/Q(i), <j 2 &3, (2 + «) 5 (1 + i) 2 ), where a 6 is a zero of the polynomial x 6 — 2x b + (Si — 51)a; 4 + 
(4i - 30)x 3 + (-2i + 755)x 2 + (—298* + 2134)a; + — 593i + 1628, is a division algebra of degree 6 and 
has a minimal discriminant. 

One of the unfortunate properties of our construction is that when we produce division algebras of a 
composite index the resulting algebras tend to have relatively large non-norm elements 7. In the following 
example we solve this problem in one specific case and show that we can always use 7 = (2 + z)(l + i). 
The method has a straightforward generalization to more common situations. 

Example 6.5: In what follows we produce the algebra A§ as a tensor product of two smaller algebras. 

Let a 2 be a zero of the polynomial x 2 + i. The algebra B 2 = (F(a 2 )/F, a 2 , (1 +i)(2 + i)) is a slightly 
modified version of the algebra A 2 of Table [TTTJ It is a division algebra with a minimal discriminant. 

The algebra £> 3 = (F(a 3 )/F, a 3 , (2 + + is a modified version of the algebra A3. Proposition 
16.41 gives us that £> 3 is still a division algebra. By considering the equation B 3 ® A3 ~ M n (F) we see 
that £>3 has the same discriminant as the algebra A3. 

Because B 2 and B3 are division algebras with minimal discriminants it follows from Proposition 15.11 
that the tensor product Aq = B3 <g> B 2 = (F(b 2 , 03)/^, a 2 a3, (2 + + i)) is a division algebra with a 
minimal discriminant. The polynomial / 6 is just simply the minimal polynomial of the generator a e of 
the field F(b 2 ,a 3 ). 

VII. An example code and some simulation results 

One of the ingredients in the construction of the perfect codes was the use of ideals in improving 
the shape of the code lattices. In [5] we did the same but for the purpose of saving energy and making 
the lattice easier to encode. We include the following simple fact (also known to E. Viterbo, private 
communication) explaining why using a principal one-sided (left or right) ideal instead of the entire order 
will not change the density of the code. 

Lemma 7.1: Let A be a maximal order in a cyclic division algebra of index n over an imaginary 
quadratic number field. Assume that the minimum determinant of the lattice A is equal to one. Let x E A 
be any non-zero element. Let p > be a real parameter chosen such that the minimum determinant of 
the lattice p(xh) is also equal to one. Then the fundamental parallelotopes of these two lattice have the 
same measure 

m(A) = m(p(xA)). 

Proof: By multiplicativity of the norm the minimum determinant of xA is equal to the absolute 
value of nr(x), so the parameter p is the unique positive root of the equation 

p n \nr(x)\ = 1. 
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TABLE III 

Conductor p of the cyclotomic field Q(C P ), 7, and the minimal polynomial /„ of the extension Q(i)(o„)/Q(i) 



n 


P 


7 


fn 


2 




(2 + t) 


x 1 + i 


3 


79 


(1 + 0(2 + 


a; 3 + x 2 - 26a: + 41 


4 




(2 + i) 


a; + % 


5 


11 


Cl + iV2 + il 


2-5 , 4 _ 4 3 _ o„2 1 , 1 

«-C J | tt' iiO | « ' - | _L 


7 


211 


(1 + 0(2 + 


x 7 + x 6 - 90x 5 + 69x 4 + 1306a; 3 + 124x 2 - 5249a; - 4663 


8 




(2 + i) 


a; 8 + i 


9 


271 


(1+0(2 + 


x 9 + x 8 - 120a; 7 - 543a; 6 + 858a; 5 + 6780a; 4 + 7217a; 3 - 2818a; 2 - 4068a; - 261 


11 


859 


(1 + 0(2 + 


x 11 + a; 10 - 390a; 9 - 653a; 8 + 52046a; 7 + 146438a; 6 - 2723930a; 5 - 1 1558015x 4 + 36326009a: 3 + 
250960565a; 2 + 385923388a; + 145865807 


13 


6163 


(1 + i)(2 + i) 


2 ,i3 + a ,i2_ 2844a; 11 - 6017a; 10 + 2908490a; 9 + 10238862a; 8 - 1340405033a; 7 -6785664624a; 6 + 
281925130086a; 5 + 1909036915713a; 4 - 21097272693753a; 3 - 192054635052100a; 2 - 
235667966495418a; + 213548387827457 


16 




(2 + 


x 16 + i 


17 


239 


(1 + 0(2 + 


x 17 + a; 16 - 112a; 15 - 47a; 14 + 3976a; 13 + 4314a; 12 - 64388a; 11 - 136247a; 10 + 422013a; 9 + 
1631073a; 8 + 41 1840a; 7 - 5840196a; 6 - 11894369a; 5 - 10635750a; 4 - 4739804a; 3 - 938485a; 2 - 
54850a; - 619 


19 


8779 


(1+0(2 + 


x w + x i8 _ 4i 58:z ,i7 + 8463a; 16 + 6281539a; 15 - 34466097a; 14 - 4291513699a; 13 + 
39454551948a; 12 + 1357034568541a; 11 - 17014625218525a; 10 - 184614267432185a; 9 + 
3035523756071878a; 8 + 10088401800577582a; 7 - 253111326110358151a; 6 - 
143208448461319868x 5 + 10612439791376560471a; 4 - 3774559232798357892a; 3 - 
220041647923912963182a; 2 + 86083932120501598139a; + 1794221202297461499641 



Let us denote this by 

p=i4ri 1/n - 

nr[x) 

On the other hand, the index [A : xA] = \N^/q(x)\ (see [20, Exercise 7, p. 131]). Moreover, [20, Theorem 
9.14, p. 119] tells us that 

\Na/q(x)\ = \N F/Q (N A/F (x))\ Rema ^ k °|iV F/Q (nr(x)")| [ ^ ]=2 \nr{x) n \ 2 = \nr{x)\ 2n . 

Hence, [A : xA] = \nr(x)\ 2n . Scaling the lattice xA by the factor p will multiply the measure of the 
fundamental parallelotope by p 2n . The claim immediately follows from these facts by calculating 

m(p(xA)) = p 2n2 m(xA) = |_L^ 2 « 2 /« u : X A] m(A) = |— ^-| 2n \nr(x)\ 2n m(A) = m(A). 

nr(x) ' nr(x) 

■ 

We remark that the same fact obviously also holds for principal left ideals of a maximal order. A way 
of using the above lemma is that we can choose the element x in such way that the left (or right) ideal 
xA is contained in the natural order. By moving the code inside the natural order we then to some extent 
recover the layered structure of the natural orders, and then, hopefully, also some of the advantages of 
the inherent orthogonality between layers. 
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TABLE IV 

Conductor?) of the cyclotomic field Q(C p ), 7, and the minimal polynomial /„ of the extension Q(V— 3)(a„)/Q(-\/— 3) 



n 


P 


7 




2 


5 


(V=3)(2) 


X 2 + X — 1 


3 




(2) 


a; 3 - 3x + 1 


4 


5 


(V=3)(2) 


x 4 + x 3 + x 2 + X + 1 


5 


101 




x 5 + x A — 40a; 3 + 93a; 2 — 21a; — 17 


7 


197 


(>/=3)(2) 


x 7 + x 6 - 84x 5 - 217a; 4 + 1348a; 3 + 3988a; 2 - 1433a; - 1163 


8 








9 




(2) 


a; 9 - 9a; 7 + 27x 5 - 30a; 3 + 9a; + 1 


11 


353 


(>/=3)(2) 


x 11 + x w - 160a; 9 - 525a; 8 + 6066a; 7 + 26034a; 6 - 48369a; 5 - 265374a; 4 - 42966a; 3 + 405001a; 2 + 
63189a; - 170569 


13 


4889 


(V=3)(2) 


a; 13 + a; 12 - 2256a; 11 + 15535a; 10 + 1555245a; 9 - 20301911a; 8 - 255557592a; 7 + 4688166666a; 6 + 
3148489502a; 5 -327998691680a; 4 +1203189132463x 3 +3781862679467a; 2 -26224493395483a;+ 
33207907136809 


16 








17 


9011 


(V=3)(2) 


2.17 + x i6 _ 4240a; 15 + 17305a; 14 + 5727403a; 13 - 41284287a; 12 - 2705219919a; 11 + 
14589308035a; 10 + 564280956214a; 9 - 1381250312443a; 8 - 51961946136288a; 7 + 
526852031838a; 6 + 1834916754576839a; 5 + 1836850197549204a; 4 - 23335163152861586a; 3 - 
34406356236297728a; 2 + 60102147038980885a; + 73569709231092527 


19 


8171 


(V=3)(2) 


2-19 + 2,18 _ 3 870x i7 + 41421a; 16 + 3724805a; 15 - 43503449a; 14 - 1437461514a; 13 + 
12225751511a; 12 + 286728047867a; 11 - 968096767438a; 10 - 28322179217822a; 9 - 
31203374649750a; 8 + 994413740064487a; 7 + 3501119135247182a; 6 - 8098862899035075a; 5 - 
59620882192114428a; 4 90513387045636018a; 3 3449524754137218a; 2 + 
73725797301678129a; + 35046894150872059 



For example in the case of the Golden+ algebra we can use the element (1 — A) 3 from the ring of 
integers Oe of the larger field E = Q(V2 + i) as a multiplier. Thus, by denoting 



M 



i-\y 

(l + A) 



we get the ideal X consisting of matrices of the form aMM x + bMM 2 + cMM 3 + dMM 4 , where the 
coefficients a, b, c, d are Gaussian integers and the matrices M, , j = 1, 2, 3, 4 are from Section HV-Al This 
ideal is a subset of the natural order Oe © uOe- 

Our code constructions are based on selecting the prescribed number of lowest energy matrices from 
a chosen additive coset of the ideal X. In order to reach a target bandwidth utilization of 4, 5 or 6 bpcu 
we thus selected 256, 1024 or 4096 matrices. In this sense we have done some coset optimization for the 
Golden+ codes, but make no claims as to having found the best coset. For the rival Golden code from 
[10] the coset corresponding to assigning all the Gaussian integers the value (1 + i)/2 stands out. This is 
because then there are 256 matrices all having the minimal energy, and more importantly because in that 
case pulse amplitude modulation (PAM) can be used to good effect. We first did some simulations using 
a PAM-type rule for larger subsets of the Golden code as well by arbitrarily selecting a suitable number 
of coefficients of the basis matrices from the set {—3/2, —1/2, 1/2,3/2} so that the desired bandwidth 
efficiency was achieved. This is a natural choice well suited for e.g. the sphere decoding algorithm. While 
we ended up having a dead even race BLER-wise at 4.0 bpcu, the Golden code lost to the Golden+ code 
by about 0.9 dB at the higher rates (see Figure [U). In the interest of a fair comparison we then tried 
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Fig. 1. Block error rates at 4, 5, and 6 bpcu. 

coset optimization for the Golden code as well. This narrowed down the gap to about 0.3 dB. However, 
the resulting subsets of the Golden code no longer have such a structure well suited to PAM. In other 
words both the rival codes must resort to the use of a code book. We have not even attempted to solve 
the problem of optimizing the code book for the purposes of minimizing BER. This also explains, why 
our performance plots only show the block error rates (i.e. the probability of decoder deciding in favor 
of a 2 x 2 matrix other than the transmitted one) rather than bit error rates. Thus, our simulations may 
also be viewed as measuring the amount of power lost, when one insists on not needing a code book. 

VIII. Concluding remarks and suggestions for further work 

We have derived a bound for the density of fully multiplexing MEMO matrix lattices resulting in codes 
with a unit minimum determinant. The bound only applies to codes gotten from the cyclic division algebras 
and their ideals. While the bound is not constructive per se, we also showed that it can be achieved for 
any number of transmit antennas, and discussed techniques leading to the construction of CDAs with 
maximal orders attaining the bound. R. Vehkalahti is preparing an even more number theoretical article, 
where these techniques are expanded. We also discussed the Ivanyos-Ronyai algorithm that is needed to 
actually find these densest possible lattices inside these CDAs, and gave as an example a construction of 
a fully multiplexing 2x2 code that outperforms the Golden code at least for some data rates. 

We have not yet exhausted the box of optimization tools on our code. E.g. the codes can be pre- and 
postmultiplied by any complex matrix of determinant one without affecting neither its density nor its 
good minimum product distance. In particular, if we use non-unitary matrix multipliers, the geometry 
of the lattice will change. While we cannot turn the lattice into a rectangular one in this manner, 
some energy savings and perhaps also shaping gains are available, but we have not solved the resulting 
optimization problem yet. Hopefully a suitably reformed version of our lattice will also allow a relatively 
easy description of the low energy matrices. This in turn would make the use of the sphere decoding 
algorithm on our lattice more attractive. 

There are also possibilities for applying these class field theoretical techniques to slightly modified 
density problems of ST-codes. E.g. it is probably relatively easy to adapt the bound of Theorem 13.21 
to the case of multi-block ST-codes. Another possibility is to study the cases, where the codes are not 
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fully multiplexing. Such situations arise naturally in an application, where the receiver may have a lower 
number of antennas, e.g. in a cellular phone downlink. 

An immediate open problem is to utilize maximal orders of the cyclic division algebra of index 2 
with center Q(u;). When looking for the example code in the previous section a natural step was to use 
LLL-algorithm for finding a relatively orthogonal basis for the lattice. That definitely aided the search 
for a good coset. In the hexagonal case this step is somewhat trickier and using a multiplier to put the 
maximal order inside the natural order only lead to a code with a disappointing performance. The best 
way of using this densest known lattice of 2 x 2-matrices is not known to us. As another open problem 
we ask, whether the discriminant bound can be broken by a MIMO lattice that does not come from a 
cyclic division algebra. We believe this to be a very difficult question. 
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